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ABSTRACT 


The radiational model used in this study computes the 
planetary albedo and absorption of solar insolation by the 
ocean's surface and by atmospheric layers for the primitive 
equation model of FNWC. Large-scale cloud parameterization 
in several layers was utilized in these computations. Solar 
insolation disposition was computed from the water-mass and 
the cloud amounts over each gridpoint. Long-wave cooling 
effects were computed using emissivity formulas after 
Sasamori at the earth's surface and over the same layers, 
and were also found to be dependent upon the cloud para- 
meterization. 

Two forms of the cloud parameterization were tested us- 
ing FNWC data over the oceanic gridpoints for 16 April 1974. 
The objective was to determine the parameterization which 
better verified the radiational balance as a function of 
latitude when compared with satellite climatology at the 
top of the atmosphere for the same general data-period in 
1969 (after Raschke, et al, 1973). The better verification 


resulted with the smaller cloud parameterization values. 
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I. INTRODUCTION 


This thesis is a study of a radiative heating parameteri- 
zation for use in the Fleet Numerical Weather Central (FNWC) 
prediction system. The study has as a primary objective, the 
analysis of the radiational and heat balance of the ocean- 
atmosphere system utilizing FNWC gridded data-fields at con- 
stant pressure levels for 16 April 1974. The gridpoint data 
were selected along four oceanic meridians, three in the 
Pacific and one in the Atlantic, with most of the gridpoints 
located in the Northern Hemisphere (Fig. 1). 

The specification of amounts of clouds in two designated 
layers has the most influence on the radiative-model disposi- 
tions (short- and long-wave). The initial specification of 
the fractional amounts of CL(1) and CL(2) are based on large- 
scale formulations developed by Smagorinsky (1960) and used 
in a similar study based on the data day 16 January 1974 by 
Spaeth (1975). 

The governing equations of the radiational transfer of 
the model have been derived by Martin (1972, 1974), who modi- 
fied both the solar and terrestrial radiation transfers to 
respond to the presence of clouds. These clouds were of 
specified amounts in the two layers, one of which is a mid- 
level cloud and the other a low-level cloud. The radiational 
model in use here has similarities to those in use in UCLA 


and NCAR General Circulation models. 
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The physical descriptiíon of the radíative model may be 
applied to any scale of analysis for which there is adequate 
resolution of the temperature and moisture data in the verti- 
cal. In the horizontal, the reliability of the data used 
here is consistent with that of the FNWC analysis to grid- 
points, and is typically reported to the nearest tenth of a 
degree with regard to temperature and dew point. The radia- 
tional computations made here are applied to FNWC gridpoints 
and are designed to make a one-hour forward-time step appli- 
cable to the FNWC primitive equation forecast model, with 
special adaptations to their o-levels, 

The Smagorinsky (1960) cloud-specification of CL(1) and 
CL(2) when used with the model of this study resulted in 
global albedo estimates that were too high as compared to 
the satellite data of Raschke (1973), thus leading to under- 
estimates of the net radiative balance at the top of the 
atmosphere and at the ocean surface. This problem suggested 
the reduced cloud-amounts el formulated by Spaeth (1975) 
and used in this study. 

The comparative results afforded by the 2/3-CL parame- 
terization gave reasonably close agreement with the radia- 
tive climatology of Raschke et al (1973), for the NIMBUS III 
period 1-15 May 1969. Oceanic and tropospheric balances 
were computed for both sets of cloud conditions (full-CL 
and 2/3-CL) for 16 April 1974. It was also necessary to 


augment the purely radiative model by including a turbulent 


1793 





boundary-layer model for sensible and latent heat transports 
over the ocean as adapted from Kaitala (1974). 

The radiation package results computed with the 2/3-CL 
cloud model compared very favorably with satellite clima- 
tology, though only one data day (16 April 1974) was used as 
representative of mid-April soundings over the ocean, both 
for the model and the Raschke climatology. However, the 
latter was based on the data period 1-15 May 1969. 

The major effect noted by the use of the 2/3-CL as con- 
trasted with full-CL parameterization was the reduction in 
global albedo and the resulting greater contribution to the 
surface net heating rates. This result is in general agree- 
Ment with recent observations by Van der Haar and Hanson 
(1969), Van der Haar and Oort (1973) and Raschke et al (1973). 
However, the major climatological discrepancies resulting 
from the specification of layered cloud amounts in the 
tropics in a manner identical to that used poleward of the 
subtropics still seems to give somewhat excessive albedo. 

It is felt that the geometry of solar radiation streams im- 
pinging upon subgrid sized cumulus cells in the tropics 
should be remodeled so as to divert a greater percentage of 
solar radiation downward, as compared with the characteris- 
tic cloud-layering and resultant reflective effects which 


seems permissible in middle and high latitudes. 
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II. DATA PREPARATION 


A. INITIAL DATA FIELDS 

The temperature and humidity data used in this study were 
arranged in the form of soundings taken along four oceanic 
meridians (Fig. 1) of the Fleet Numerical Weather Central 
(FNWC) Northern Hemisphere analyses on 16 April 1974. Compu- 
tation of radiational dispositions and of other heat budget 
terms are made at these gridpoints. Oceanic locations for 
these computations were chosen because: 

l. The constant 0 surfaces (where О = < ) of the FNWC 
primitive equation system are close to being constant pres- 
sure levels. 

2. The maritime-area heating rate computations are like- 
ly to be representative of the month of April 1974 as con- 
trasted with computations made for 16 April 1974 over a 
corresponding set of gridpoints over land. 

The three meridians (and their respective number of sound- 
ings) selected over the Pacific Ocean were located at 125W 
(25 soundings), 170W (25 soundings) and 145E (17 soundings). 
The Atlantic Ocean meridian was 35W (26 soundings). This 
method of selecting "soundings" along the indicated meridians 
of the FNWC polar stereographic map, made it unnecessary to 
employ spatial interpolation between original data gridpoints 


along the meridians. Data along line 3 in the Pacific was 
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FNWC polar stereographic grid and meridians (lines 
1, 2, 3, and 4) selected for study. The longitudes Á are 


shown for each meridian as well as the extent considered of 
each meridian. 
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not extended southward of gridpoint (9,55) because they fell 
over land masses (New Guinea and Northern Australia) where 
the surface temperatures and other sounding features were 
unrepresentative of the oceanic values. 

The gridpoint soundings in the form of Table I(a) were 
taken from the original FNWC 63-by-63 surface analysis, as 
well as of the vertical distribution of T(P) and of dew 
point depression at five other standard pressure levels up 
to and including 400 millibars. These original soundings 
were then modeled into a radiative sounding (Table I(b)) cor- 
responding to FNWC primitive equation prediction levels up 
to 100 mb (see Fig. 2). 

Data from 0000GMT, 16 April 1974 were used for the Paci- 
fic soundings while 1200GMT, 16 April 1974 data were used 
for the Atlantic soundings. These times were used so that 
the analysis times would correspond as closely as possible 
to local solar noon in each of the indicated areas, 

At each gridpoint selected, the original humidity data 
was given in the form of five dew point depressions for the 
analysis levels from 925mb to 400mb. Since the standard in- 


strument level vapor pressure (e ) was missing from the 


alr 
original FNWC analysis at all gridpoints, e irp was approxi- 
mated using the FNWC field of e , a computed value of vapor 
pressure at about 20 meters above mean sea level (MSL) in 
the turbulent boundary layer. FNWC values of T. and e, over 


the ocean had been made available using an operational 


planetary boundary layer model detailed by Kaitala (1974). 
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Table I(a). Example of typical FNWC sounding for gridpoint 
(1,1). The humidity parameters between 925,.., 

400mb are dew point depressions, but at the surface and top 

of the constant flux layer (9999%*) are vapor pressures in mb. 


Pressure (mb) TWE) Humidity Parameters 
Surface, 1000 mb 25 776 SS 27.4 mb 
925 18,3 T 2 “G 
850 14.6 1.2020 
700 56 2:4 6 
500 -11.0 205 C 
400 -22.7 - 1,4 C 
300 -37.6 
250 -46.5 
200 -56.3 
150 -67.3 
100 -80.6 
9999% 5 23.5 


* 
Code 9999 indicates data taken from the top of the constant 
flux level (CFL) of the FNWC initial data progranm. 
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Table I(b). Example of the corresponding radiative sounding 

with mixing ratio listed at odd k-levels (Fig. 2). 
Additionally, water vapor and CO, absorber masses, cloud 
amounts, CL(1) and CL(2) (after Smagorinsky (1960)) as these 
parameters are modeled in the radiative theory. 


155. Absorber Masses Smagorinsky 
sure T(°C) PATIO water er CO, : Cloud Amounts 
(mb) g/kg (gm/cm ) (cm/em ) CL(1) CL(2) 

1000 25.60 17.10 
900 16.67 12.56 1.000, CL, 
800 11.84 2435 45.53 
700 5.60 6.43 
600 =L; 57 3.36 893.053 
500 -11.00 27763 0.900, CL, 
400 -22.70 1.41 34069 TE. 
300 -37.60 0.68 
200 -56.30 5775 b... 
100 -80.60 0.04 2575 36.0 / 
0 -80.60 . 00 3.75 143.39 
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исиге 2, Five-layer radiative sounding used in this study. 
Levels are identified by their values on the k-scale, while 
layers are identified by their level boundary indicies in 
parentheses, e.g. (8,10). Pressure-scaled water vapor and CO, 
mass increments AU and 4C, respectively are integrated rel- 
ative to the surface and the resulting U and C are carried in 
the model at even levels. The temperature T is retained at all 
levels. Amounts of clouds CL(1) and CL(2) in the layers shown 
have been parameterized for consideration of their radiative 
effects. 
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To perform radiative calculations, it is necessary to 
have water vapor and CO, absorber masses and cloud amounts, 
CL(1) and CL(2) at certain required k-level boundaries 
(Fig. 2). All soundings in this study start at sea level 
with the approximation of surface pressure m = 1000mb. 
Therefore, the k-levels correspond approximately to the 
FNWC levels Pi IA “00. «ses, 200., 100., 0.0mb 


corresponding respectively to O, Zu E09 2,0, 1,0080: 


B. INTERPOLATIVE PROCESSING TO k-LEVELS IN RADIATIVE 
SOUNDINGS 


1. Temperature Profiles 
The gridpoint temperatures were listed at each manda- 
tory level of Table I(a) between 1000., ..., 100mb. The tem- 
perature was assumed to be isothermal from 100mb to p=0.0mb. 


The temperature T was set equal to the FNWC listed sea-sur- 


10 
face temperature. The radiative sounding temperatures for 
the remaining k-levels were obtained from either their corres- 
ponding listed temperature-level or by a three-point Lagran- 
gian interpolation scheme (Eq. 2-1, Spaeth, 1975), to level 
k when the listed FNWC temperature-profile did not include 
the value TE 
2. Moisture Profiles 

The moisture from the original FNWC soundings (Table 

I(a)) were converted into mixing ratios at each of the 


original sounding levels. The near-surface vapor pressure 


(taken as 55 was used to calculate the mixing ratio at k=10 


2.47 





as described by Spaeth (1975): 


E21 377 Ce /2:000) OE) 


The remaining ratios q were calculated from the original 
FNWC sounding (Table I(a)) dew point depression data after 
Spaeth (1975, Eq. 2-3). The computed q-values were subse- 
quently interpolated to k-levels using the previousiy men- 
tioned three-point Lagrangian procedure applied to successive 
q-values in the original sounding. Resulting values of the 
interpolated q-values are shown for the case of the radia- 
tive sounding at gridpoint (1,1) (Table I(b)). 

Due to the fact that most radiosonde humidity data 
for p < 300mb is either unreliable or not available, a power- 


law extrapolation formula (Spaeth, 1975) 


tk ск) 
ds 500 (2-2) 
was used to obtain q-values at k=3,2,1. 
This extrapolative procedure for obtaining q-values 


at high levels was first suggested by Smith (1966), utilizing 


the least squares estimation of A according to 


6 A oo 
ии: 2 
il 5 7 >; 
A ee R "7 (2-3) 
Mae, yx x. 
21 i 


q. 
with y. = log — and x. = log(P./P.). The Smith method was 
J. qs 1 hə 
then further tested by Spaeth (1975) for the vertical scale 
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O MEMA NVC mi tial data soundings., The resulting correla- 
tion coefficient 5 vas generally found to lie in the range 
+95 to .99. Such high correlations indicate that the pro- 
cedure described by Spaeth (1975) for determing upper atmos- 
phere q-values has realistic estimation value for the 
vertical scale involved in the radiative sounding. 

3. Pressure-Scaled Absorber Masses 

The pressure-scaled water vapor absorber mass u in a 
layer (Fig. 2) was calculated for the five odd numbered k- 
levels by using the computed mixing ratio values (Eq. 2-7, 
Spaeth, (1975)). The equation for computing the integrated 
pressure-scaled water vapor mass is given by the algorithm 
described by Eq. 2-8 of Spaeth. 

A similar algorithm was used in the computation of 
the carbon-dioxide scaled absorber mass (Spaeth, 1975, Eq. 
2-10,11). To be noted here is the fact that CO, absorber 
masses have been stated in terms of N.T.P. "pressure-scaled" 
cms, since the CO, absorption coefficients are generally 


2 


stated in terms of this CO, mass unit. 


C. CLOUD PARAMETERIZATION 

The relative humidities, and thus the saturation vapor 
pressure at levels k=5 and k=9, are used in the calculations 
of the fractional cloud cover CL, and CL, in layers (4,6) 


and (8,9) respectively, as depicted in Fig. 2. The fractional 
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cloud amounts for these two layers were parameterized using 


the following equations (after Smagorinsky, 1960): 


CLE(DR=>2.0 (RH(5)) – 0.7 (2-4(a)) 


CLE2)M= 3.35(RH(9)) – 2.0 (2-4 (b)) 


Smagorinsky's parameterization of CL(1) and CL(2) permits 
cloudliness fractions of 1.30 and 1.333 respectively with 
RH=1.0. Values of CL greater than 1.0 were considered 
initially by Smagorinsky (1960) to suggest the amount of 
supersaturation which accompanies precipitation, 

All operational radiative models in use at the present 
time limit CL < 1.0. In the radiative study developed here 
no opportunity is afforded to deduce precipitation-rates 
and/or supersaturation amounts. Hence, the Smagorinsky for- 
mulations of CL, and CL, were reduced in each case by the 
multiplicative factor 2/3 so that neither fractional cloud 


cover could exceed unitv. Thus alternative cloud-cover 


fractions 


' 
CL 1 


' 
CL 2 


273 (CL,) (2-5(a)) 


2/3 (CL,) (2-5(b)) 


vere considered for purposes of estimating the shift in the 
global radiative balance as a result of altering the cloud 
computations in the radiation model, 

220717 reduction inherent in CL',> CL", relative to 


Smagorinsky"s CL CL vercd initiated in this study in an 


1? 2 


30 





attempt to tune cloud amounts to 57 radiation results in 
closer agreement with recent satellite climatology of 
Raschke et al, (1973). This satellite climatology suggests 
the use of smaller cloud amounts than that specified by 
Smagorinsky. Fractional cloud amounts were considered to 
be functions of large scale effects only. Therefore, small 
scale convective activity, seasonal conditions (except for 
the soundings) and latitudinal effects were not considered 
in specifying the reduction factor of Eq. 2-5. The reduced 
cloud-parameterization was introduced here for estimating 


large scale radiational effects only. 


D. CLOUD-AREA COVERAGES 

Since Eqs. 2-4a,b or 2-5a,b gave the fractional-coverage 
of the gridpoint area by the appropriate cloud-type, the 
gridpoint area may be thought of as broken into random frac- 


tional segments of size 
W(0,0) .= (1-CL, ) (1-CL,) — (2-6a) 


wherein there is a combination of clear-over-clear segments 
in the layers. Similarly, the gridpoint area has the frac- 


tional area of cloud coverage 


ти си (2-6b) 


of upper-cloud amount CL, overlying lower-cloud amount CL 


2 28 


Likevise the combinations of cloudy over clear and clear 


1 





over cloudy , by layers, may be visualized as occurring with 


the respective weights 


Il 


W(1,0) CL, * (1-CL,) (2-6c) 


1. 


ап4 


W(0,1) 


(1-CL,) кет, (2-64) 


2 


Regardless of whether the full-cloud amounts CL; > CL, of 
Smagorinsky (Eq. 2-4), or the 2/3-CL amounts of Eq. 2-5 were 
utilized, it was useful to carry, for radiation computations 
of each sounding, the relative weights or fractions of the 
gridpoint area exposed to the specified cloud-layer combina- 
tions. Henceforth, the symbols denoted by W(0,0), W(1,1), 
W(1,0), W(0,1) and as given by Eqs. 2-6(a,b,c,d), suggest 
overcast (1) or clear (0) cloud amounts in the indicated 
layers (Fig. 2), the first index 1 or O referring to layer 
CL,» and the second to CL, 

The usefulness of this computational device will be 
clarified in Sections III and IV, where the procedures for 
the terrestrial and solar radiational computations are dis- 


cussed and the results are summarized over the set of 


soundings. 
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III. TERRESTRIAL RADIATION 


E THEORETICAL AND EMPIRICAL BASIS 

Empirical formulas were developed by Sasamori (1968) for 
flux emissivities in the atmosphere associated with computa- 
tions for the heat balance requirements of the NCAR General 
Circulation Model, Sasamori derived the empirical emissivity 


formulas for water vapor and CO, by comparison with the 


2 
theoretical values built into the Yamamoto Radiation Chart 
(1952). This chart has proved to be quite accurate for com- 
putational checking of the Sasamori emissivities and was used 
by Spaeth (1975) and Warner (1974) as a systematic guide for 
integration of the radiative transfer formulas developed by 
Martin (1972, 1974), who adapted the Sasamori formulas to 

the particular layers of interest in the gridpoint computa- 
tions of the FNWC primitive equation model (Fig. 2). 


The essential long-wave flux formulas required for use in 


the FNWC heating package are the following: 


Fio“ net IR flux at earth, k=10 

Fo m net IR flux at level k=6 

F,* = net IR flux at level k=2 

F610 = net IR flux divergence in the layer (6,10) 


F26 = net IR flux divergence in the layer (2,6) 
To get the IR flux divergences in the layers (6,10) and 


(2,6), the differences FL *-F. q * and F.-F.” must be computed. 
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The detailed scheme for making such computations with various 
combinations of cloud cover CL(1) and CL(2) is similar to the 
technique developed by Martin (1974) and as reproduced in 
detail by Spaeth (1975). 

In order to make IR net-flux calculations along the path 
of integration, there must be a physically sound representa- 
tion of the emissivity (£ ve as a function of both vater 


vapor and CO, absorber masses in layers along the sounding. 


2 
For a further discussion of the emissivity formulas used in 


the quadrature scheme, refer to Spaeth's Appendix A (1975). 


B. NET FLUX FORMULATION 
1. At Level k=10 

The radiative sounding as depicted in Table I(b), was 
computed as the combination of parameters U(k,10), C(k,10) 
and Ti for each required level, k=10,8,...,2,1,0. Cloud para- 
meters CL(1) and CL(2) are also listed at each gridpoint and 
in general are non-zero. The grid area was then considered 
to be composed of areal E Lione (weights) defined in 
Eqs. 2-6(a,...,d) and denoted by the symbols W(0,0), W(1,1), 
Peele, 0), W(0,1). 


The composite net flux F (CL) > CL.) at level k=10 


ln 


at each gridpoint is then constructed by using the appropriate 


weight factor to multiply the reference net flux F__* compu- 


10 


tations defined for the special cloud-cover cases 


* * * * 
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so that 


* = < * 
F,¿*(CL, ,CL,)=59(0,0)F, ¿*(0,0)+W(1,0)F, ¿*(1,0) 


+HW(0,1)F,,*(0,1)+9(1,1)F,,* (1,1) (3-1) 


Spaeth (1975) has listed these reference net flux formula- 
tions in his Eqs. 3-6, 7,8. Using the definitions of W(0,0), 
1770), V(0,1), V(1,1), Fig* (CL, CL.) can be shown to assume 


the form 


Fjo* (CL, »€L,) = [1-CL (2)11(B,p7B¿)-»5l€, (8,10) (B,-Bg) 


+ (E, (8,10)4 E, (6,10)) (B2-B,)]) 


+ 


2 ... O SUE, (6,10) 
+ E (4,10)) (B,-B,)+(e£. (4,10) (3-2) 
+ £ (2,10))(B,-B„)+(eE (2,10) 


m 1 0 BB 777 E 


+ 


CL(2){(B,)-By) [1-.5€, (9. 10) | Q<: 


Here 


B, = 1.170403 x 107? ə (323) 


is the Stefan-Boltzmann blackbody flux in langlies per day. 


Further, € 10) is the combined water-vapor 


215 


and CO, emissivity along the path from level 10 to level k. 
This emissivity is considered by Sasamori to be temperature 


independent for T > 210K, whereas 7: represents the 
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temperature dependent emissivity r R LOr TE TOK 
(see pp. 136-137, Spaeth, 1975). 

The reference net fluxes Ејо" of Eq. 3-1 are asso- 
ciated with (1) clear skies in both layers, (2) overcast in 
the upper layer only, (3) overcast in the lower layer only 
and (4) overcast inboth layers respectively. 


2. Ne tölbülux Fa 


6 


The formula for FGACCL,, 


Spaeth (Eq. 3-10, 1975) in a manner analogous to the deriva- 


CL.) has been developed by 


tion of the weighted Fig” The result is reproduced after 


Spaeth (1975, see pp. 43). 


X = zın la 1777060 (BB ve (4,6)(В.-В,) 
FC Le (Sə 50 FE  (2,6)XB,-B,)+(£. (2,6) 
A (Beeb tee. (((0,6),7,)*B, 1} 
+(1-CL (1)) (1-CL(2))((B, ,7-Bg)[1.-.5(e, (6,8) 
ETERO (3-4) 
EI E Ba)" 


IR, (6,0)! ј 


3. Net Flux F,* 


The net flux at level k=2 can be calculated in a 
similar fashion to Fio and Few, with the formulation of the 


individual reference net fluxes F,* as previously described. 
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The formula for F,*(CL,,CL,) is then developed analogous to 


* * 
that for Fig (CL,,CL,) and Fg (CL, ,CL,). 


F *= [1-CL(1)](B2-.5[e (2,4) (B,-B,)+(e, (2,4) 


+ €v (2,6))(BL-B,)+(€ ` (2,6)+€ (2,8))Х 
(BgrBg)t+E (1,2)(B,-B,)+È  ((0,2),T,)4B.]) (3-5) 
+(I-CL(1)) (1-CL(2))“((B, -B,)[1-.5(£. (2,8) 


+ 2771 7 000 0 516  (2,5)(B,-B.) 


+ 7777707 7” 0077775) È 


A new parameter, the total outgoing long-wave radia- 
tion to space (FF2) can be readily defined from the expres- 
sion for Fo” by setting to zero the terms of Eq, 3-5 
representing the dovnvard flux through level k=2. These 
terms are iust those denoted by the symbols Zo 
and eve (1929 (B2-8)9. If the long-vave radiative model, in- 
cluding its cloud parameterization, is realistic, the modeled 
outgoing flux to space (FF2) should compare reasonably with 
that observed by NIMBUS III satellite during the same period, 
after Raschke et al (1973). The most nearly comparable 
satellite climatological period was 1-15 May 1969. These 
satellite observations were obtained from the NIMBUS III 
satellite atlas of Raschke et al (1973) for the May dates 


indicated and used for comparison with the model values of 


FF2, point for point (at 5°latitude intervals) along the 
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samel four meridians of Fig. l. Also the latitudinally- 
distributed mean satellite fluxes were compared with the 
corresponding averaged values of FF2. The results showing 
the comparisons of the model computations of FF2 with the 
corresponding satellite observations are listed by latitude 


in Table V, 


C. APPLICATIONS TO HEAT BALANCE COMPUTATIONS 
l. At Upper and Lower Boundaries 
In order to compute the heat balance at the top of 
the earth-tropospheric system the composite F,* and the 
total insolation absorbed below level k=2 at 7 gridpoint 
are required. Moreover a radiative balance is computed at 
the surface for each gridpoint using the earth's absorbed 
insolation (Section IV) and the composite net flux loss Fio*: 
In order to consider the possibility of a heat 
balance at the earth's surface, rather than merely a radia- 
tive balance, it is necessary to add a term representing com- 
bined sensible and latent heat transfer across the air-sea 
Interface (Section V); This latter transfer rate at each 
gridpoint was adapted from the FNWC primitive equation model, 
after the discussion of Kaitala (1974) and based upon the 
FNWC data of 16 April 1974. 
2. Intermediate Levels 


Consideration Of the radiation balance in the atmos- 


phere requires computation of the long-wave cooling effects 
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caused by the flux divergences F26 and F610. F26 and F610 


are defined as 


= kou * a 
F26 F, F Ç (3-6) 


= di = * = 
F610 Fe Fio (3-7) 


where the symbols (2,6) and (6,10) indicate the layer 
boundaries involved in the balance considerations for the 
indicated layer, 

In order to compute F26 and F610 at each gridpoint, 
complete listings of IR net fluxes Fo > Ес̧”, Fo, associated 
with each set of reference cloud amounts (0,0), (1,0), (0,1), 
(1,1) have been computed at each gridpoint together with the 
weighted set of fluxes F, (CL »CL,),F¿*(CL, ,CL,), 
F,*(CL,,CL,). These cloud-weighted values have been con- 
structed using the weighting scheme of Eqs. 2-6, 3-1, 3-4 and 
3-5. A sample gridpoint printout of the net fluxes has been 
included in Table II. 

The determination of the heat balance by layers has 


been deferred to Section VI, where the results are displayed 


in cross-sectional form for the 2/3-CL case. 


D. STATISTICAL RESULTS AND COMPARISONS 
1, Net Flux F. ,*(0,0) 


A statistical test of the F__*(0,0) numerical results 


10 


was conducted by utilizing a linear regression program for 


relating the predictand F,,*(0,0) against the simultaneous 


10 


values of X"Bio and Be as predictors, using all 93 
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gridpoint data as samples. The regression program was based 
on the BMDOZR in the Biomedical set of programs (Dixon, 1973). 
The result was cast in the form of the well-known Brunt net 


flux equation (1932): 


EH ON bye) (3-9) 


where 


4 
— de 


is the Stefan-Boltzmann blackbody flux at the surface tempera- 


ture Tio and e is the surface vapor pressure in mb. The 


multiple correlation coefficient ти was higher 

than 0.99 in the 16 April case, as was also true for the data 
periods of Spaeth (16 January 1974), and Beahan (16 July 1974), 
(Table III). 


Since Fio* may be defined as 


* = = = 
Fio Bio Fa (3-10) 


Eq. 3-10 can then be solved for Fi 


UB hub , (3 17) 


for the clear-sky case. It must be recalled that the compu- 


tation of Fi included both the effects of water-vapor and of 


CO Since there is no direct relationship of Ve upon the 


2 
quasi-constant absorber mass of CO, , thewradiative  éffect/oftf 


the latter constituent in Eq. 3-10 must be included in the 


term involving (l-a) Bio: 
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A cross-seasonal comparison of the coefficients l-a 
and -b of Eq. 3-11 is presented in Tabie III, comparing the 


clear=sky results of F, for 16 April 1974 with those of Spaeth 


d 
and Beahan at intervals of three months earlier and later, 
respectively. 

The high multiple correlation coefficients indicate 
that an oceanic version of the Brunt downward flux equation 
is valid for each tested season, Рә depends most strongly 


upon B while only slightly on B,pYe: This is an apparent 


10 
result of the nearly constant relative humidity over the 
large oceanic regimes. Moreover, the weak dependence of Ej 
upon B,pYe is most evident at the data-time 16 July 1974, 
when the coefficient -b is smallest. This characteristic is 


due primarily to the relatively weak gradient of Ve in the 


radiative soundings aloft. 


16 January 1974 16 April 1974 16 July 1974 


l-a ә .6557 .8129 

-b .0224 .0209 -.0094 
Бу. .4500 25577 51 
Ri .9955 .9969 .9950 


Table III. Cross-seasonal comparisons of the coefficients 
l-a and -b from Eq. 3-11, The mean values F 
: : d 
from the Brunt equation and the resulting multiple 
correlation coefficients have been included. 


2 e ladification of Fio” for Cloudiness CL, CL 


A second statistical test was performed based on a 


2 


possible relationship between the cloudy-sky and the clear-sky 
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< i = * k 
cases of Flo 5 05 ra o Y Fio (CL), CL,)/F j o (0,0) was 


used as the predictand in this case. The predictor CL was 


the total opaque cloud cover at each gridpoint, defined as 


(after Quinn, 1971 and Spaeth, 1975): 
Cin ac) — CL€22:(1-CL(1)). (3-12) 


This regression vas developed utilizing Flo” (CL,,CL,) in 

the numerator, computed with the amounts CL(1) and CL(2) as 
calculated from Smagorinsky's formulation, and then utiliz- 
ing Fig*(ch',,cL',) computed from the 2/3-CL parameterization. 


The general form of the regression formula required is shown 


in Eq. 3-13 below, with CL given by Eq. 3-12 above. 
* = * É 2 
Ep” (CL) x F. 70,0) .11-d(CL)1 . (3-13) 


A cross-seasonal comparison of the coefficients d and other 
SIH S tics for the 2/3-CL and fulMCL casesM using the results 
of this study and those of Spaeth and Beahan, is presented 
in Table IV. There seems to be little clear-cut variation in 
the statistics of Table IV between the results of adjacent 
seasons, and little statistical preference between the full- 
Gitana the 2/3-CL cases. 

S 0 correlation fecertrerenetsioi Table IV for 
all seasons considered, indicate the general capability of 
the IR radiative model employed in this study to account for 


the effective net radiation at the surface. 
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J6 January 1974 16 April 1974 16 July 1974 


2/3-CL Full-CL 2/3-CL Full-CL 2/3-CL Full-CL 
d ./5 . 75 «73 “75 . 79 .01 


Ry ©9936 . 9916 . 9945 . 9942 — 0 . 9916 


Table IV. Cross-seasonal comparison of the coefficient d 

from Eq. 3-13 and the resulting multiple corre- 
lation coefficients for both the 2/3-CL and the full-CL 
cases. 


The mean values of Fig” over the 93 gridpoint sound- 


ings for 16 April 1974, as applicable in Tables III and IV 


were as follows: 


Full-CL Statistics 2/73ECL Statistics 
= — Zij — = 
* = ï * = . 
Flo (CL) .0639 ly min Fig (CL) .0868 ly min 
ONO) = .1515 ly mim | MANN 005) nin: 
10 10 
a 7733 CL = .5888 


The superior bar symbol denotes the sample mean. 
These results indicate that the surface net flux is decreased 
by the respective ratios .578 with CL = 277( and by ¿27 
with CL = .5888. These results indicate that the downward 
f bux Ber) is 5/.8 percent larger at the surface than 
F¿(0,0) when full-CL amounts are used in the cloud parameteri- 
zation model. In the 2/3-CL parameterization the mean down- 
ward flux at the surface BALDI is increased by only 42.7 


percent when compared to F,(0,0). 


AEC COMPARISONS OF FF2 WITH SATELLITE CLIMATOLOGY 
Comparison was made of computed-model values of FF2 with 


satellite measurements of total long-wave flux to space 
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(after Raschke et al, 1973) for the NIMBUS III period of 

l-15 May 1969. EEZ is used here instead of F,* since the 
downward flux at k=2 is not measured by the satellite. The 
computation of FF2 was previously discussed. FF2 for 16 
April 1974 was computed for both full-cloud and 2/3-CL cases 
at all 93 gridpoints. These FF2 values were then inter- 
polated to whole multiples of 5 degrees of latitude along 
each of the four meridians and the results were then averaged 
to get a mean latitudinal distribution value for each 5 
degree latitude interval considered in this study. 

Table V shows both of the latitudinally distributed FF2 
cloud-model values compared with those extracted from 
Raschke et al (1973). Raschke's results were obtained by 
interpolation from charts using the same oceanic meridians 
as those used in this study. Latitudinally distributed total 
Opaque cloud cover CL by Eq. 3-12 are also listed for both 
the full and 2/3-CL cases at each latitude. In the bottom 
line of each column in Table V is listed the cosine weighted 
Mean of each set of the column values for the listed latitude 
range (see Eq. 7-2 for the equation defining "weighted 
Average"). 

It should be noted that the Raschke results do not corres- 
pond to specifically known CL values. Both cloud cases give 
FF2 results which are reasonably close to the values reported 
by Raschke, especially the means and between 5N-25N. Neither 
the full-CL nor the 2/3-CL model can be conclusively selected 


as preferable, based on Table V results alone. The limitations 
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of the comparisons made here are obvious, when it is re- 
called that between latitudes 20S-5S and between 60N-65N 
there are fewer than four meridional lines available for 
computing the listed zonal values in Table V. For all other 
latitudinal average values, four meridional lines were used 


in the averaging. 


F. COMPARISONS OF CROSS-SEASONAL RESULTS OF Flo” 


Table VI depicts the latitudinally-distributed values of 
Fio* and of CL at the earth's surface obtained using the 2/3- 


CL parameterization. A seasonal comparison of the model com- 


puted F,.* is made with data obtained from Spaeth (1975) for 


10 
the 16 January 1974 and from Beahan (1975) for 16 July 1974 


cases. It is clearly shown that F_.*(CL') is a decreasing 


10 
Function of cloud cover. There is a clear-cut tendency in 
each season for a maximum value of Fio” to be located in the 
subtropics (latitudes 15N-25N). Also there is evidence of a 
high latitude (55N-60N) minimum Fio” associated with a con- 
Sr ration of maximum cloud cover CL', The most outstanding 
variation is the transition in the Southern Hemisphere lati- 
tudes (20S-10S), which has large cloud cover in both the 
January-April monsoon period but only small CL' in the July 
data period. The same relative conclusions resulted when 
the full-CL parameterization was examined for cross-season 


comparisons, but the obtained results are not included in 


Mable VI, 
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ma. FF2 F2(RAS) Cloud-Amounts 
Chul l=—CL) (2/3-CL) (Raschke) Full-CL 2/3-CL 


205 sı e205 2 .41 .985 . 842 
15 . 2688 «2962 2 . 983 828 
10 .. 3005 270: . 41 . 944 275 

23220 . 3383 . 40 . 842 .647 

25216 . 3381 77 . 878 205 

ə 3499 5012 535 . 829 2.1 
10 43711 . 3801 37 5072 . 639 
15 .3714 . 3804 7 5 .548 
20 25: .3815 39 . 788 . 536 
25 5 766 2517 7 . 642 .473 
30 . 3463 35 39 5 2262 .446 
35 . 3290 ә3373 . 3 6 .634 .520 
40 . 3099 . 3206 33 .630 .468 
45 . 2914 571 . 34 «630 .467 
50 232463 .2893 . 34 2577 .674 
DIO . 2665 2793 37 eg . 555 
60 £5: 22530 231 2 962 . 816 
65N 51615 52065 0 255 57: 
Wt avg 3264 206 . 3740 796 . 35956 


Table V. Comparison of zonally averaged longwave flux to 
space, FF2, as found by this study for both full- 

Manda 2/3-=CL cloud cases for 16 April 1974, and F2(RAS) 

by Raschke et al (1973) based upon NIMBUS III measurements. 

Also included are composite cloud-amount fractions used in 

the two cloud models of this study. Flux values in ly min-1, 
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16 January 1974 16 April 1974 16 July 1974 


Lat Fjg*($ -CL) $ -CL F,g*(È -cL) É -cL F) ($ -CL) 2 -CL 
205 078 2759 074 . 842 . 151 . 269 
15 ə 079 ./25 «067 «8328 6135 . 306 
10 . 086 .634 «071 ə /52 «135 . 258 
5 . 090 . 536 . 081 .647 .093 698 
0 . 062 .556 2073 . 659 . 100 665 
5 «086 57: . 080 591 . 0915 .534 
10 .10/ . 376 . 086 . 639 7177 .460 
15 115 . 342 . 092 .548 . 103 .430 
20 . 120 .279 . 085 536 .094 .431 
25 . 120 3:20 .094 .473 . 098 337.5 
30 elil әз ЕТ .098 .446 092 .425 
35 2177 382 “17 . 520 . 092 260 
& 0 . 108 257. .116 .468 «+ 085 . 3 3:7 
45 . 067 . 740 . 100 .467 ¿071 . 566 
50 «069 279 . 080 ə 074 . 064 2475 
55 . 068 . 749 . 108 .555 .038 70 19 
60 . 106 .538 .054 . 816 .085 2550 
65N «095 0.0 .082 . 796 061 5: 
Wt Avg .096 0514 . 086 .596 .092 .405 
Table VI. Cross-seasonal comparison of the zonally-averaged 


aate lues at the surface Es, and cloud coverage 
based on 2/3-CL parameterization for the six-month interval 
25776 january 1974 to 16 July 1974, Flux values in ly mə mm... 
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IV. SOLAR RADIATION 


A. PARTITION OF SOLAR INSOLATION 

The solar constant assumed in this study at level k=0 
(top of atmosphere) was 2.00 hy nəə (Joseph, 1971). This 
flux was depleted by 4% to account for the attenuation 
caused by oxygen and ozone above the tropopause. This left 
the value S = 1,92 Ey uu. at level k=2 to be used in this 
Study as the effective solar constant. 

Eq. (4-1) was then utilized to compute the effective 


solar insolation at the tropopause (k=2) as follows: 


mez) = 5 (E)? Cos z (4-1) 
m 
where S = effective solar constant at k=2 
Cos z = cosine of the zenith angle for the Julian date 
used 
r/r, = ratio of the actual earth-sun distance to the 


mean earth-sun distance for the Julian date used 
in this study. 


The Smithsonian Meteorological Tables gives the ratio r/r. 
and the solar declination ô for 16 April, 0000GMT, as these 
values are employed in Eq. 4-1. These values are, respec- 


tively: 


li 


= 1.00333 
T> 
m 


6 374833 degrees of lat. 
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6 is used in evaluating the cosine of the solar zenith 


angle, given by 


Cos z = Sin $ Sin Š + Cos Q Cos ô Cos h (4-2) 


where Q is the latitude and h is the hour angle of the sun 
for the meridional data lines. For example, Fig. l makes 


it clear that for 


line l h = 55° 
line 2 h = 10° 
late ses A h = 35° 


at the times of the synoptic charts (0000GMT and 1200GMT, 
as applicable). Sin 6 was computed according to the standard 
polar stereographic projection formula applicable to the FNWC 


base chart given by 


re = [(1-32)° + (3-32)°} 
Sin Q aren L 9 ee (4-3) 
Tr + {(I-32) + (3-32) ] 
where a = 973.752. Thus Sin Q assumes the following func- 


tional form in terms of the ENWC grid-coordinate I (Fig. 1): 


Lines 1,3,4 
2 
ünü. = nn. 215 S 
973.752+2(32-1) 
Line 2? 
77 72(32-1 


2 
557 5 (4-4(b)) 
973.752+(32-I) 
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or, conversely I is given in terms of Q by 


Bes 8). 05 Y] (4-4(c)) 


J q Sin Ó 


505 0 — Ò (4-4(d)) 


1 + Sin Q 


Lines 1,3,4 I 


Line 2 I 


25577 “5 for Lines 1,2 
IF 0 ToT” Lİn€e 3 


IN. Line 4 


For lines 1,2, and 3 the gridpoint soundings correspond to 
0000GMT 16 April when solar noon occurs at the 180th meri- 
dian. Line 4 gridpoint soundings correspond to 1200GMT, 

16 April, when solar noon was at the Greenwich meridian, 

A very simple partition of solar insolation was utilized 
ers study after Joseph (1971). It consisted of dividing 
the insolation F(2) into two parts at level k=2. One part 
was considered to include all wavelengths A > .9 um where 
absorption by water vapor and carbon dioxide bands are the 
most prevalent attenuation processes in clear air. This part 
of the solar spectrum was termed the F(A) energy and con- 
sidered subject to water-vapor absorption but not to Rayleigh 
scattering. For those wavelengths À < .9 Um, absorption of 
the solar insolation energy by water vapor was considered 
negligible. This part of the solar insolation was denoted 
7 Se pestivwe of the fact that it was subject only to Ray- 


leigh scattering attenuation in clear air. The two partitions 
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are formulated after Joseph (1971) as follows: 


F(A) 5 497 F (2) (4-5(a)) 


F (S) “051 E (2), 4 (4-5(b)) 


In this study, the introduction of two cloud decks produced 
cloud-reflectivity effects upon both the F(A) and F(S) 

solar energy insolations. However, in the clear areas around 
any gridpoint only the absorption-attenuation applies to the 
F(A) insolation, while only Rayleigh scattering-attenuation 


applies to the F(S) insolation. 


B. DISPOSITION OF F(S) INSOLATION 

Im the disposition of the F(S) insolation, Joseph (1971) 
determined that Rayleigh scattering reflectance to space by 
clear skies (after Coulson, 1959) could be effectively approxi- 


mated by least squares in the following form 


G (R) 


27755 2250 RN [log ə sə) (E 
O 


where P_ = Io mb. In Eq. 4-6, MP. = 1 in view of the 
fact that mean sea level pressure T is close to 1000 mb, 


Also 
-1 
Sec z = (Cos z) 


with Cos z given by Eq. 4-2. 
e ace albedo O(G) is another reflective parameter 
utilized in this study. Over oceanic areas the following 


бот о (Сб) after Gates et al (1971), was utilized: 


ac) = max {.06, .06 + .54 (.7 - Cos z)}. (20) 
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L Clear Sky Case 


In the clear sky (0,0) case the F(S) insolation was 
subjected to both Rayleigh scattering reflectance a(R) and 
to surface reflectance a(G). Considering the likelihood of 
a succession of multiple reflections between earth and atmos- 
phere, the F(S) insolation actually penetrating the earth's 


surface after scattering is given by 


1810(0,0) = F(S)[l-a(R)][l+a(R)a(G)+...(a(R)a(G))” 


+ 


...]*(1-a(G)) (4-8(a)) 


that is, by 


IS10(0,0) Ми 711-616) 1/ 11-0(R)(G)1 09), 


2 Cloudy-Sky Cases 


In the three cases in which clouds were present, F(S) 
insolation absorbed by the ground at each gridpoint was com- 


puted using the following equation (after Arakawa, 1972): 


JS100C1;1) 


F(S) (1-R(1)) (1-R (2)) (1-a(G)) 
* 1-[R(1)R(2)+R(2)a(G)+R(1)a(G) 


Bee. (48:9) 


As indicated by the notation (1,1),denoting CU OSO 
Equation 4=9 is the formula used in calculating F(S) insola- 
tion absorbed by the ground in the case where overcast clouds 
Би ирси Dot levels of Fig. 2. Also in Eq. 4-9, con- 
stant cloud-reflectance values were chosen, namely R(1) = .54 


for the mid~level clouds between k=4 and 6, and R(2) = .66 
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for the low-level clouds between k=8 and 9, Both cloud- 
reflectance values are as suggested by C., D. Rodgers (1967). 
For all the other cloud cases, the following changes 
were applied to Eq. 4-9. In the (1,0) case (CL(1) = 1.0, 
CL(2) = 0.0), the desired result is obtained by setting 


1570770 in (4729), from which it follows that 
IS10(1,0) = F(S)(I-R(1I))(1-a(G))/[1-R(1)a(G)]. (4-10) 


In the case (0,1), one sets R(1) = 0.0 in (4-9) so that (4-9) 


simplifies to 
IS10(0,1) = F(S) (1-R(2)) (1-a(G))/[1-R(2)a(G)]. (4-11) 


Note that with a cloud overcast present, the Rayleigh clear- 
sky scattering a(R) does not appear in Eqs. 4-9, 4-10 or 4-11, 
weis included in the cloud reflectances R(1) and/or R(2). 
3. Composite F(S) Insolation 

Eqs. 4-8, 4-9, 4-10, and 4-11 were utilized in the 
computation of the cloud-weighted F(S) insolation penetrating 
the earth's surface considering the areal-weights of the cloud 
donmbrnations denoted by (0,0), (1,1), (1,0) and (0,1) about 
a gridpoint. The resultant F(S) insolation penetrating the 


earth's surface denoted by IS10 is therefore expressible as 


FS PO CJ (1) CL(2)) IS10(0,0) W(0,0) 


5855101 1) V€1,I) 
+ IS10(1,0) W(1,0) (4-12) 


SINO AO NEO. 1). 
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Here the weighting factors W(0,0), W(1,1), W(1,0) and W(0,1) 
are computed in Eqs. 2-6a, 2-6b, 2-6c, and 2-6d respectively. 
Note finally that the part of F(S) insolation reflected to 
space is found by subtracting IS10(CL(1), CL(2)) from F(S). 
Table VII lists the results of radiative sounding as it 
appears for the disposition of F(S)-insolation. The indivi- 
dual computations of 1510 аз they apply for the possible 
overcast-clear layer cases are made under the heading "IS10." 


The difference 
REFS = F(S) – 1910 (4-13) 


in each case represents F(S)-insolation reflected to space 


while 


REFS 


STRAN = OS 


(4-14) 
has been computed as that portion of the F(S)-insolation 
incident at the sea surface just prior to transmission by 
the surface. Note that no absorption in air has been in- 
cluded in the computations of Table VII, and that the only 
absorption permitted is that implicit in IS10. Finally at 
the bottom of each column, e.g., IS10, the composite value 
has been computed by means of the weighting scheme of 


Eiq . 4-12, 
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CL(1),CL(2) Weighting IS10 REFS STRAN 


0,0 0.0 22000 751 5216 
ino 0.0 .2539 .3517 .3076 
0,1 .1013 5190217 .:135 .2327 
1,1 .8997 1340 .4656 .1696 


Composite F(S) values „1452 .4603 .1760 


Table VII. A sample listing of values of F(S) insolation 
(ly min-1l) computed at gridpoint (1,1) using 
equations detailed in Sec, IV, 
C DISPOSITION OF F(A) INSOLATION 
The fractional portion of the solar insolation subject 
to absorption by atmospheric water-vapor and carbon dioxide 
are covered in the following subsections. 
l. Clear-Sky Case (0,0 
The Manabe-Moller absorptivity function provided the 
necessary absorptivity values for the key layers in this 


case. The form of this absorptivity function is 


Memo = eek) Sec z]*"93 Cais 


Here a is the absorptivity applied to the pressure-scaled 
water vapor mass between levels 2 and k (Fig. 2) along the 
zenith slant-path angle z. The resultant absorbed insola- 
tional energy in the particular layer (2-6) is then given 


by the Manabe-Moller relation 


MOM 0005011060 6) see z)1”77” . (4-16) 
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The two layers of interest in which absorption was computed 
were (2,6) and (2,10). The absorbed insolation in the layer 


(6,10) was then computed by 
OE AZ LO) - A(2;6) >» (5217) 


Water-vapor mass above level 2 was assumed negligible in the 
F(A) disposition of the solar insolation. 

By subtracting A(2,10) from F(A), the direct trans- 
Mission of F(A) insolation impinging at the earth's surface 
was determined. The transmission of F(A) insolation is then 
further reduced by the transmissivity (l-a(G)), after sur- 


face-reflectance which leads to the earth-absorbed result 


20016): 


IA10(0,0) = F(A)(1-.271110(2,10)Sec zl) 
(4-18) 


The transmitted energy impinging upon the earth iust prior 


to absorption is 
TRANA(0,0) = IA10(0,0)/[1-a(G)] . (4-19) 


Za Cloudy Cases 


In order to compute meaningful dispositions of F(A) 
insolation in cloudy-sky cases, cloud reflectivities and 
25: 0:5 iviZies after €, D, Rodgers (196/) were utilized. 
The reflectance-values used here are different from those 
suggested by Rodgers for the F(S) wavelengths. The cloud 


reflectivities used here are RA(1) = .46 and RA(2) = .50. 
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In this case there are also cloud-absorptivities to be con- 
sidered. These were taken as A(1) = .20 and A(2) = ,30, 
respectively. In the following discussions the cloud con- 
ditions are considered totally overcast in each of layers 
indicated when the notations ((1,1), (0,1) and (1,0)) are 
utilized. 

Schematic representation of the computations per- 
formed in the vartous cases ((1,1), (0,1), (1,0)) are dis- 
played in Fig. 3. This figure indicates the theoretical 
parameters required in computing the disposition of incoming 
PCA)-insolation from level k=2 to the earth's surface (k=10), 
The equations which relate to the parameters in Fig. 3 and 
the other cloud configurations are listed in Appendix B of 
Spaeth s study (1975). 

The insolations, A24, A46, A68, A89 and A910, etc., 
represent the contributions to the insolation absorbed in 
the layers (2,6) and (6,10) involved. Symbols F2, F4, F6, 
F8 and F9, etc., depict the streams of insolation passing 
through the indicated level. A vertical arrow implies the 
direction if insolation passage, i.e., + denotes downward 
insolation, t upward-reflected insolation, and yy downward- 
reflected insolation. The absorption quantity A(6,8)++, for 
example, indicates absorbed energy remaining in (6,8) from 
a downward reflected beam. 

Since multiple reflections occur between the earth's 


surface and a cloud base, or between cloud-layers in this 
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F(A) 
F2? 
2 
A24? A241 
4 An. RA(1 | F41 
A(1) , A46 
ó / İ RA(1) | 
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1A10 


Figure 3, schematic representation of F(A) insolation 


disposition in the case of two overcast layersi(c a s e(11)), 
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model, the insolational amount reaching the lowermost reflect- 
ing surface is affected by the downward reflected beams, 
Additional insolational amounts remaining after more than 

two reflections are negligible, and are not included in the 
model computations. Thus, note that insolation reflected up- 
ward from a lower interface, either cloud-top or ground, to 
the base of a higher cloud deck was not subjected to absorp- 
tion by or transmission through the cloud. This simplifica- 
tion resulted in slightly reduced F(A) insolation-reflectance 
еомерасе:, 

A final point to note is that there are three con- 
tributions to the absorption of solar insolation within a 
layer between any two reflecting surfaces. In Fig. 3 these 
contributions are identified by arrows which indicate the 
portion of path being crossed. 

For each of the possible overcast-cloud combinations, 
the transmitted F(A) insolation arriving at the earth's 
surface may be defined using the following notation (TRANA) 


as 


TRANA(1,1) = F10+ + F10++ (4-20) 
TRANA(1,0) = F10¥ + F10¥4 (4-21) 
TRANA(0,1) = F10+ + F10++ (4-22) 


where the right side parameters of (4-20), (4-21), and (4-22) 


are derived in Appendix B of Spaeth (1975). 
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The F(A) insolation Absorbed by the earth in each 


cloud-case considered is derived from 


MAOKI DS F104 (12a(G)) + F104+ 4-23) 
Podio FIO, Cl-a(G)) + F10++ (4-24) 
wooly = Plot (l=a(G)) + F10++ (4225) 


In the last three equations, the quantity F10¥+¥ is small 
enough in each case, that no further reflections from the 


earth were considered. 


3% Composite F(A) Layer-Absorptions and Surface- 
Absorption Insolation 


As has been previously discussed, the standard grid- 
area weighting scheme of this study was applied to obtain 
composite values of the absorbed F(A)-insolation in key 
layers and also within the earth's surface. The weighting 
factors applied to the corresponding overcast-combination 


absorption quantities provided the following composite results: 


A26(CL(1), CL(2)) = A26(0,0)W(0,0)+A26(1,0)W(1,0) 


+ 


A OO ит) (4-26) 


+ A26(1,1)W(1,1) 


A610(CL(1), CL(2)) =A610(0,0)W(0,0)+A610(1,0)W(1,0) 
+A610(0,1)W(0,1) (4202) 


+A610(1,1)W(1,1) 


1501000007 GL(2)) =IA10(0y0)W(0,0)+IA10(1,0)W(1,0) 
+IA10(0,1)W(0,1) (4-28) 


+IAIO(1,1)W(1,1) 
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The weighting factors W(0,0),...,W(1,1) were first listed in 
Eqs. 2-6(a,b,c,d). Also A26(0,0), A610(0,0), IA10(0,0) are 
KA a ni ear sky case (0,0) about each gridpoint by 
Eqs. 4-16, 4-17 and 4-18, respectively. 

1: results tor the absorption in layers (2,6), (6, 
10) and at the surface are shown in the following table for 
бә i ource ioudicasesmr'(0;.0); (1,0), (0,1), (1,1)] and for the 


composite cases based on the weighting factors (also listed). 


CL(1),CL(2) Weighting A26 A610 AI10 REFA TRANA 


0. 0 0.0 520 000065 OS04.:1348 .1751,7069 
1950 0.0 SD 00306 20482 „3517 23774 
0,1 1043 25:21 1499770238 4140 . 2697 
deg . 8987 M2970717: 0104 (4656 .1880 
Composite values 2177770777 .0118 .4603 .1962 


Table VIII. A sample listing of gridpoint values (ly (n 
of F(A) insolation computed at gridpoint (1,1) 
using Eqs. 4-26, 4-27 and 4-28. 
In the computational scheme indicated by the entries 


of Table VIII, the reflected F(A) insolation to space has 


been depicted by the symbol REFA, and its values follow from 


REFA = F(A) - A26 - A610 - IALO (4-29) 


Whereas the TRANA dispositions are given by Eqs. 4-19, ..., 
4-22, respectively, or by its weighted-mean value in the 


case Of TRANA-composite. 
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4. Absorptivit ABA) by Layers 
Here the (fractional) absorptivity as well as the 
actual insolation values absorbed in the layers are con- 
sidered., In the computation of absorptivity, which is frac- 
tional absorption, the total undepleted insolation at the 
top (k=0) is used as a base. The following equation was 


caia zed in this calculation: 
FADJ = 2.00(r/r 77 Woe 2 E (4-30) 


The absorptivity of the troposphere ABA was computed from 
the ratio of the insolation absorbed in the troposphere to 
the insolation incident at the top of the atmosphere rather 


than at k=2: 


ABA = 2500 + AC6710 775 


FADJ ° 


D. ALBEDO (ALB) OF THE EARTH-TROPOSPHERE SYSTEM 

In considering the planetary albedo, the reflected inso- 
lations of the earth-troposphere system in both the F(A) and 
F(S) insolational regions must be recalled by the program. 
Thus REF is computed at each gridpoint as the sum of the re- 
flected insolation energy in F(A), denoted REFA in Eq. 4-29 


2417716 re 1ected part of F(S) (denoted REFS in Eq. 4-13): 


REF = REFS + REFA . (4-32) 


Finally the planetary albedo is related to FADJ through 


REF 
2 Даи 
— PADS ° 
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E. COMPOSITE ABSORPTIVITY (ABG) BY THE EARTH-SURFACE; 
COMPOSITE ATMOSPHERIC TRANSMISSIVITY (ATRAN) 


l. Absorptivit ABG) of Earth 
By summing the weighted values of F(S) and F(A) por- 
tions of the incoming insolation entering the earth, the 
total insolation absorbed at the earth's surface was com- 
puted. This quantity when divided by the extraterrestrial 
insolation gave the fractional absorptivity (ABG) of the 


earth's surface. The equation for ABG was 


IA10 + IS10 
5256 = si (4-34) 


where IA10, IS10, and FADJ were defined previously by Eqs. 
4-28, 4-12, and 4-30, respectively. 
2. Transmissivit ATRAN) of the Troposphere 


Also computed was the total insolational energy TRAN, 


e 


incident at the earth's surface just before absorption by 


the surface. This calculation is given by 
TRAN = TRANA + STRAN . (4235) 


Here STRAN = [ISIO/(I-a(G)] was previously defined in Eq. 
4-14, and a(G) was given in Eq. 4-7. TRANA has also been 

215£.7 0705: the weighted value of TRANA(0,0), TRANA (1,1), 
TRANA (1,0) and TRANA (0,1) given by Eqs. 4-19, 4-20, 4-21, 
and 4-22. Note also in justification of STRAN that the four 
cases for IS10 of (4-8), (4-9), (4-10), and (4-11) each have 


the common factor (1-a(G)) in the numerator and therefore 
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each transmitted F(S) insolation component available at the 
earth just before absorption needs only be divided by 
(1-a(G)). TRAN may thus be viewed as the total insolational 
energy incident at a pyrheliometer located at earth. The 
(fractional) transmissivity of the troposphere (ATRAN) is 


then computed from 
ATRAN = TRAN/ FADJ . (4-36) 


Note finally that the major dispositions of the total in- 
solation at the indicated map times have now been identified 
by the fractional values, ALB, ABA or ABG, and ATRAN, repre- 
senting the reflectivity (albedo), absorptivity of air or 
earth, or atmospheric transmissivity as the case may be. 
3. Computational Check 
The computational scheme utilized in this model was 

checked continuously by summing the fractional values ALB, 
ABA, and ABG at each gridpoint. The value in each case must 
total .96, since as previously noted the attenuation of 
solar insolation was taken as four percent as it passed 


through the stratosphere. 


F. STATISTICAL ANALYSIS 

Im order to substantiate some of the computations per- 
ned in this section, several of the most important itens 
computed were statistically analyzed using linear regres- 
sion computer programs from the BMD set of statistical pro- 


grams (Dixon, 1973). 
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lee Gila Sky Cases 
Uu 6 0.0,2), vhich is the clear sky case of 
albedo, as the predictand, and logio Sec z as the predictor, 


the following best-fit equation resulted 


EE UTUĞZ) 07734 + .4007210g, Sec VA 


RA” . 9624 (4-37) 


where Ru represents the multiple correlation coefficient. 
This result was as expected over the ocean where 


a(G) and a(R) involved logarithmic dependence on Sec z. The 


sample average values of ALB and 108419 Sec z were 
ALB(0,0) = .1358 log, Sec ZAS O -a 


Other clear-sky regression tests made use of the water-vapor 


mass path length (M) defined by 
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Hü (U Sec 5 log,,(U Secaz) 4 (4-38) 


This parameterization of water vapor mass is similar to that 
developed by Hanson (1971), who used a similar M as a pre- 
Its empirical formulations of ABA for both clear 
and partiy cloudy sky cases. In this study the two regres- 
In attempted using M as a predictor were with ABA(0,0,M) 
and ATRAN(0,0,M) as predictands. The best-fit equations 


resulting vere 


ABAKO TOM) = 2.09572 +. .03136M (4-39) 
Ry = «9875 
ATRAN(0,0,M) = .79520 - .03182M (4-40) 
557 .9488 . 
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Smesabsorptivity in the form ид. by (4-39) agrees closely 
in form with that obtained by Spaeth based on 16 January 
1974 data and with that of Quinn based on pyrheliometric 
daşa (1971), 

The means of ABA(0,0,M) and ATRAN(0,0,M) for the 


original 93 gridpoint sample vere 


ABA(0,0,M) = .13459 


ATRAN(0,0,M) = .75576 


M 1.23945 (gm cm” 2) 142 ; 


дә. Statistical Relationships Between ALB, ABA and 
ATRAN in the Cloudy and Clear Sky Cases 


In addition to computing the clear-sky value of ALB, 
it was also possible to compute the composite cloud value of 
ALB. A regression was formed showing the relationship be- 
tween the ratio ALB(CL(1),CL(2))/ALB(0,0) as the predictand 
and the total opaque cloud cover (CL) as the predictor. CL 
in this model is specified by Eq. 3-12. CL, it was felt, 
gave a good e ana aa of the effective cloud-cover by 
the two layers of cloud amounts CL, and CL, by Eq. 2-4 used 
in this model. As explained earlier in this study, a one- 
Ae dieron of both CL(1) and CL(2) after initial determi- 
nation by Eq. 2-4 was also tried in all computations involving 
cloud amounts. This test represented an attempt at tuning 
the cloud model for radiative calculations, and gave results 


in closer agreement with the latest literature on satellite 


reflectances. Thus all subsequent regressions of ALB, etc. 
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upon CL, after Eq. 2-4, will also include those best-fit 
equations formulated using a CL' based upon Eqs. 2-5 and 
3-12. -A small letter subscript 'a' in the equation number 
will indicate those computations using 2/3-CL amounts while 
a small letter subscript 'b' will indicate those utilizing 
full-CL cloud model amounts. 

The first regression equations tested were those for 
ALB(CL) and ABA(CL) for the respective cases of 2/3-CL and 
full-CL cases. The results are given in (4-41(a)), (4-41(b)). 
The symbol RM once again signifies the multiple correlation 


of the statistical regressions, 


ALB_(CL(1)CL(2)) = ALB(0,0)[1 + 4.7154CL 
cani da) 


Ey = .9640 


ALB, (CL(1)CL(2)) 


oo 7 CL 
50770: (4-41(b)) 


Ru = „9670 


The sample means of the values in Eqs. 4-4la and 4-41b are 


Eq. 4-4la Eq. 4-41b 

ZE 01(2)) - .3624 ANEL (an), Ch(2)) = 4473 
ALB(0,0) = 1256 ALB(0,0) = .1358 
GL = .5876 CT = .7746 


The high value of the correlation coefficient in Eq. 4-41 


shows strong dependence of ALB upon CL in either cloud 
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paorumererization case, Also note that ALB(CL(1),CL(2)) 
values computed here are clearly higher than reported by 
Raschke et al, (1973), where ALB = .25 for essentially the 
same set of gridpoints during the NIMBUS III period 1-15 
May 1969. The 2/3-CL cloud parameterization gives albedo 
results which are closer to Raschke's results, suggesting 
that the full-CL amounts of Smagorinsky (Eq. 2-4) are too 
high for radiational computations. 

A second regression was developed between ABA(CL) 
/ABA(O,0) versus CL. The following best-fit equations 


results: 


Il 


ABA (CL(1),CL(2)) = ABA(0,0)[1 + .5873CL 


- .1018CL”) (4-42(a)) 
Ry = «8956 
ABA, (CL(1),CL(2)) 


ABA(0,0)[1 + .6526CL 
- .1481CL”) 60 


Ку = „9012 


Mean statistics for this regression case are 


Eq. 4-42a Eq. 4-42b 


ABA (CL(1),CL(2)) 51776 ABA, (CL(1),CL(2) = ,1894 


.1346 


.1346 ABA(0,0) 


ABA(0,0) 


CL “1250 


x (5876 CL, 


From the means of (4-42), it is seen that the model specifies 
an increase in atmospheric (solar) absorptivity with increas- 


177 0000 cover. This result is in agreement with studies 
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of Warner (1974) and Spaeth (1975) and with the earlier 
study of London (1957). 

An analogous statistical regression was then de- 
veloped for the cloudy-sky transmissivity relative to the 


clear-sky transmissivity. The resulting regressions were 


ATRAN (CL(I),CL(2)) = ATRAN(0,0){1 - .6581CL] 


Ry =" 00902 , (4-43(a)) 


ATRAN, (CL(1),CL(2)) ATRAN(0,0)[1 - .6833CL] 


Ry = .9986 (4-43 (b) ) 


with sample mean values as follows: 


Eq. 4-43a Eq. 4-43b 


ATRAN (CL(1),CL(2)) 


«4653 — ATRAN (CL(1),CL(2))= .3589 


ATRAN(0,0) ,7558 ATRAN(0,0) = .7558 


The parameter ATRAN(CL')/ATRAN(0,0) given by (4-43a) is in 
very good agreement with the statistical result of Savino- 
Angstrom from earth-based pyrheliometry at Canton Island, 


which was 


Q = од (E J655CL) 


27 CU7::“:(7:::6)- The coefficient of CL in this result is 

known to be a slowly decreasing function of latitude. 
550015 the parameters ALB(CL')/ALB(0,0) and 

ABA(CL')/ABA(0,0) behave consistently with respect to increas- 


1077707 wan more realistic mean values for ALB(CL') and 
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ABA(CL') than for the corresponding full-CL parameterization. 
Hence, in the mean, the statistical analyses performed seem 
to show agreement with observational results of other investi- 
gators. Comparative computations of the atmospheric trans- 
missivity, particularly in the 2/3-CL case appear to give 
better agreement with ground-based pyrheliometry. 

A major test still to be made is the comparison of 
ehe Patitudinal-distribution of ALB (CL) and ALB (CL) with 


that of the satellite albedo of Raschke et al (1973). 


G. ALBEDO COMPARISONS WITH PUBLISHED RESULTS 

The tropospheric albedo computations of the solar-inso- 
lation model of this section have been presented for both 
the full-cloud and tvo-thirds CL cases, respectively. These 
albedos have been interpolated to whole multiples of 5” lati- 
tude between 20°S and 65°N, These computations were made 
for each of the four oceanic meridians and the resulting 
albedos averaged across the meridians are presented as a 
function of latitude in Table IX for both the full-CL and 
two-thirds CL cases, where comparison is also made with 
satellite albedos of Raschke et al (1973), for the NIMBUS 
III period 1-15 May 1969. For sake of consistency the 
Raschke albedo data have been extracted from atlas presenta- 
tions to the 5 degree latitude grids, on the same four 
meridians, before applying the same averaging process. 

The model results under the headings "FULL-CL" and "2/3- 
CL" correspond to the composite cloud amounts previously 


listed in Table V according to the Eq. 3-12. 
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The most unusual results to be observed in Table IX are 
the small albedo-values (* .2) between 20S to 20N, after 
Raschke,.as contrasted with the corresponding values (in the 
range .3 to .47) by the 2/3-CL model, and more markedly with 
those by the full-CL model (values in the range .43 to .59). 
The cosine-weighted mean albedos (see Eq. 7-2) are as 


follows: 


FULL-CL, ALB = .46 
2/3-CL, ALB = .37 


RASCHKE-ALB = .25 


At every latitude ($ < 60) under comparison the 2/3-CL 
albedo computation is closer to Raschke's observations than 
is the full-CL value so that there is no question in accept- 
ing the 2/3-CL over the full-CL case. At latitude Q = 65, 
there is evidence of added surface-ice reflections contribut- 
ing to the Raschke value of planetary albedo. For example 
at latitude 65N Raschke's results indicate ALB = .5, whereas 
the 2/3-CL model computations with a rather high cloud cover 
MA. yan is 47. Part of this discrepancy with the 

solar radiation model at high latitudes lies in the non- 
inclusion of a surface-ice reflectance parameter. 

To summarize, apart from the region covered by sea-ice, 
m5? 7 a marked preference for the 2/3-CL over the full- 
CL model in giving albedos closer to those observed by 
1:70 76£ the same general period (1-15 May 1969) as 


reported by Raschke et al (1973) and Von der Haar and Hanson 
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ALBEDO MODEL 
Cloud-Amount 


Lat. — (Full-CL)(2/3-CL) (Raschke) (Full-CL)(2/3-CL) 
205 .568 .473 “07 "985 . 842 
15 559 .462 617 .983 .828 
10 .512 .418 .18 «944 2752 
5 .446 .363 “20 , 842 . 647 
0 2475 .374 a] .878 . 659 
5 443 . 347 22 .829 . 591 
10 .438 .337 20 .823 .639 
15 .434 .333 7: . 799 .548 
20 .430 2331 .18 .788 .536 
25 . 387 .314 .20 . 642 Bun 
30 57 . 298 .24 . 582 446 
35 .361 .316 .29 .634 .520 
40 155 .311 . 28 .630 .468 
45 414 .336 .30 .630 .467 
50 . 499 .415 .33 «877 . 674 
55 443 .383 .40 . 719 .555 
60 .562 .493 77 .962 .816 
65N .609 ‚474 .50 .995 .796 
Wt Avg .458 .373 . 24 .796 .596 


Table IX. Comparison of planetary albedo as found by this 

study for both full- and two-thirds cloud models 
for 16 April 1974, and by Raschke et al (1973) based upon 
NIMBUS III measurements. Also included are the globally- 
weighted mean values at the bottom of each tabular column 
and the composite cloud amounts for both cloud cases. 
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(1969) and others. The primary difficulty even with the 
2/3-model occurs in tropical and subtropical latitudes 

where Von der Haar and Hanson (1969) have discussed the 
reality of even smaller than suspected cloud-covers and 
smaller values of the resultant global-albedo. At these 
latitudes (tropical and subtropical) it seems doubtful that 
a large-scale CL-parameterization of the Smagorinsky type 
will be applicable. Rather some method of specifying the 
distribution at gridpoints of isolated cumulus towers may 

be appropriate. Finally, it is open to question whether 
vertically-structured convective cloud elements will have 
the high reflective capability attributed to clouds which 
are depicted as existing in horizontal layers with a con- 
stant reflectivity .5 < R < .65 as specified in the radia- 
tion model. It may well be that the vertical shafts between 
convective elements may be a more efficient focus for solar 
energy to the earth than is perceived by ascribing to the 
cloud amount a constant cloud albedo value, A more realistic 
approach may be to reduce the two-cloud layering in the 
tropics to clouds existing iín a single layer. 

At any rate with the two layered cloud-model considered 
in this study, the only feasible choice is the 2/3-CL para- 
meterization. This selection is carried over into the com- 
putations of Section VI. A further test of the 
cloud-parameterization effect on the computation of net 
radiation of the earth-atmosphere system, by comparison with 


that observed by satellite, is made in Section VII.D. 
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V. SENSIBLE AND LATENT HEAT TRANSPORT 
AT THE SEA-AIR INTERFACE 

A. GENERAL PURPOSE 

The model used here to describe the turbulent transports 
of sensible and latent heat across the sea-air interface was 
basically that already in use in the operational FNWC primi- 
tive equation model. This model has been discussed in detail 
by Kesel and Winninghoff (1972) and by Kaitala (1974). The 
primary purpose of the adaptation of the turbulent flux model 
in this study was to test the magnitudes and directions of 
these transfer rates in comparison with those of the radia- 
tional transfer model discussed here in Sections III and IV. 
The final purpose of the inclusion of the turbulent flux 
model was to conduct a diagnostic computation of the heat 
balance within the oceanic and atmospheric layers for 16 April 
1974 under the combined operation of the radiational and 


turbulent flux models, 


55002 NDS PEED COMPUTATION IN THE TURBULENT FLUX MODEL 
At all of the 93 gridpoints tested, the geostrophic wind- 


speed Vg was computed by use of the standard formula 


AZ 


_ AZ _ -1 
ME nu 2 (d/m) (5-1) 
with 
Е Av : 
ren 2 (523) 
_ 1.86603 
1551077 
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and 


AZ = y (2,290 + (2,-2,)° (5-3) 


In Eq. 5-1, d = 381.0km is the nominal FNWC grid-spacing 
(true at 60°N) and (22) is the true earth spacing used in 
the centered-difference computations of the contour gradient 
75:57 Fig. 3). 

In Eqs. 5-1 and 5-3 AZ is the contour-difference, cen- 
tered in most cases on the gridpoint under consideration. 
Hovever, at the map edge gridpoints of the meridians, İt vas 
only possible to compute a forvard-differenced version of Vg 
7775717 4). Thus at edge-point "0", Eq. 5-4 below was 
used. The spacing 2 in Eq. 5-4 is then taken as one-half 


that in (5-1), so that 


v7 aa; (5-4) 
with 
22 KOUN + 520 (5 25) 


All geostrophic windspeeds were then converted to surface 
windspeeds using the following empiricism due to Langlois and 


Kwok, (1969) 
an = 8 Vg f 2.2 (5-6) 


where VE = surface windspeed (mps) and Vg is the 1000mb geos- 
trophic windspeed. The factor 2.2 ín (5-6) is an empiricism 


introduced to include the effects of gustiness upon the 
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measured (steady-state) windspeed. 
In the interval of + 10 degrees of latitude, the value 
of the Coriolis parameter was arbitrarily set at f =+.25 x 


FB... following Kaitala (1974), while for |ọ| > 10 de- 


10 
grees the actual value of f for the particular latitude was 


used. 


C. SENSIBLE HEAT TRANSPORT 

The details of the adaptation of the FNWC model for sen- 
Sible heat transport Hr utilized in this study follow the 
description of Spaeth (1975). The formulation of Hr is 


given below in terms of bulk parameters of the surface and 


of the planetary boundary layer (10,8) of Figs. 2 and 5: 


Hr = Pip En Cy V va — T.) 
(5-7) 

руд = CASA 

m0 cd 0 
vhere 

С̧ = .239 cal 050060 
and 

Су = 1.4 х TO Cara Weiler and Burling, 


1067) 


The equation utilized in this study to compute T. may be 


shown to be (Spaeth, 1975) 


2, 5 x 
C VT + —у— - 2ER W 
= — (5-8) 
by ~ 2 K 
AZ x uns 


ad 





where 





1000. 286 
to ano (5-9) 

Sa — (5-10) 

l+a E) 

9 “10 

R 
= = 4 10 

DE 29 210 e T (1n g) (5-11) 


In Eq. 5-10, the turbulent parameters K* and a* have been 


determined empirically by FNWC 
K* = lio u 


5x 10” ee) > ; 


ax 


Do - de in Eq. 5-7 is the temperature change for the thin 
layer between the surface and the top of the constant flux 
layer (cf., Fig. 5(a)). The formula for T. deduced here re- 
quired a constant value of sensible heat Hr then decreasing 
linearly with pressure to Hp = 0 at k=8, and remaining zero 
above k=8. This linear decrease of Hp with pressure is com- 
patible with the condition of constant convergence of the 
sensible heat flux per gram in the turbulent boundary layer 
(x,8). Thus, the constant value of the convergence of sen- 
sible heat as given by Eq. 5-7 within the layer (8,10) may 
then be considered to be applicable at level k=9 as depicted 
Wi Mo), after Spaeth (1975). 

After computing T. by Eq. 5-8, Hr was then computed at 


dp ont using the bulk transfer equation 5-7. The 
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Figure 4. Section of FNWC polar stereographic meridian il- 
lir ing the method of obtaining contour gradients in the 
Vicinity of a gridpoint. The values of contour-heights 
correspond to the 1000 mb D-values in the neighborhood of 
the gridpoint. 
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gridpoint values of Hr are used later in calculations of 
atmospheric and oceanic heat balances in Sections VI and 


VIII. 


D. EVAPORATIVE HEAT TRANSPORT 

Similar to the case of sensible heat convergence, the 
latent heat transport by turbulence was assumed to be sub- 
ject to constant evaporative flux convergence throughout 
the layer 800-1000mb. The total amount of latent heat re- 
moved by evaporation from the ocean surface was essentially 
as modeled by Kaitala (1974) using bulk-transfer theory. It 
should be recognized however, that the те ләт and con- 
sequent realized latent heat may occur at arbitrary levels, 
i.e., in general at levels higher than k=8, Thus, in fol- 
lowing the FNWC model, we have not introduced E as a heating 
rate at level k=9, 

The basis of the turbulent latent heat model adopted is 


shown in Fig. 5(c), where 


E (900) 1/2 E (Bulk) (5-12) 


and where the surface layer bulk transfer of latent heat is 
given by 


E = S E(Bulk) = Lp o Cp V. (q, - q) : (5-13) 
Evaporation at level k=9 may also be computed using 
2557570 


eu) — g) Gals? 


10 


where K is given in Eq. 5-10. 
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In (5-13) and (5-14) L is the latent heat of vaporization. 


in SOUL sab. Im (5-13) and (5-14) qo and q, are the mix- 
ing ratios at the levels k=9 and at x, respectively. The 
factor 1/2 which appears in the right side of (5-12) re- 
emphasizes the concept of Bon Sen latent-heat convergence 
ieee bevel (x,8), (cf., Fig. 5(c)). 

The solution for ud is obtained by combining Eqs. 5-12, 


5-13 and 5-14 in the form 
Ом —- (5-15) 


where AZ = 24 - 210" 

Applications of Eqs. 5-13 and 5-15 make possible the 
computation of E at each gridpoint. As in the case for Тә 
ene factor of 2K oí (53-]5) was táäken to conform to the value 
K as used by Kaitala (1974) so that the turbulent transport 
results obtained here will be compatible with those computed 
operationally by the FNWC primitive equation model. As 
already noted, this was accomplished in this study by using 
a value of 2K* = 2 x OS SS as equivalent to the K*- 
value of Kaitala (1974). Thus computations of E by Eqs. 
5-13 and 5-15 were equivalent to that using Kaitala's scheme 


(1974), the only change being that the vertical distribution 


2215770136 in Fig. 5(b) is slightly different from that 
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of the constancy of E everywhere in the planetary boundary 


layer, as postulated by Kaitala (1974). 


E. TURBULENT HEAT TRANSPORTS OVER AN ICE-COVERED OCEAN 
Though this study did not have any data taken from an 
ice-covered ocean, it should be noted that Spaeth (1975) de- 

vised an acceptable method of handling this realistic con- 
sideration, which would be applicable if our computations 
extended beyond the poleward limit of this study which was 
established at $ = 65N. 
F. LARGE-SCALE TURBULENT HEAT FLUX ACROSS THE SEA-AIR 

INTERFACE 

As shown in Section VI, the combination of E + Hr turns 
out to be a heat loss mechanism at the ocean surface, as 
computed at each ae... The values of Hr which repre- 
sent the sensible heating rates for the atmospheric layer 
(8,10) were computed for each gridpoint along the four 
meridians. The values of Нр in general turned out to be 
small and negative in low latitudes, reflecting the computed 
negative values of T - T. assalsssreported by Spaeth (1975). 


š 


Im the mid-latitudes, where . - T. > 0, Hp was positive and 
therefore contributes positively to the net atmospheric tem- 
perature changes in the layer (8,10). The latent heat term 
(E) was dominantly large and positive relative to Hr through- 
out this study. Thus, regardless of possible inconsistencies 


in the sign of H combined values of E + Hr seemed to give 


E” 


reasonable heat-loss rates at the ocean surface at all lati- 


pudes mas also found by Spaeth (1975). 


82 





"Tone, STE 


+queqsSuoo əuq fo doz ayi FO TaAet ƏU SF ,X, afOoN *(9) (4461) ETEITEY 19332 хптј 


jeou-JualeT pue *(q) Apnys STU} UT pasn tleou-quener “(e) Apnas Stu UT posn Iesu 


STATSUaS JO souUazZiIJAUOJ fo WOTANATJIISTP au Jo uorqequəsədadəqa otieuouoe 


"G əqnsıq 
(9) (4) (E) 
— (3) TAN 1524 (sə) jo 40142044 ——— (3) xnj|j joey juayoj jo 402024 A — (JH) "PI 4034 0/q17u35 jO 101/3043 
01 G' 0 ol or ç o ol ot с 0 








a2u206s2au05 ON 02090519400 ON 














10314 изәјој 
jo 


әзмәбзәлиог ON 


may justo] 


jo 


10044 əlqisuəş 
jo 


ozuobsoAuoj 


oo} juejoy e2vabsanvay 





jo 


әзмәбзәлио-у 


83 





In the large-scale sense E + Hp is included as a positive 
contribution to the heat budget of the tropospheric column, 
even though it is not known precisely where the latent heat 
source E will be realized. The use of the equations of con- 
servation of mass, momentum and of water-vapor mass, would 
resolve the latter question, but this was considered beyond 
the scope of the present work. 

The turbulent flux computations for Hp and E outlined in 
the foregoing subsection V,(C). and (D) represent minor modifi- 
cations of the FNWC turbulent boundary layer model which has 
been appended to our radiation model (for the 2/3-CL case) 
for heat balance computational consistency tests, that is, 
to determine whether the net radiative flux at the surface is 


of the proper magnitude, relative to E + Hn. 
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vr, MERIDIONAL CROSS-SECTIONAL DEPICTION 
OF THE HEAT-BALANCE COMPUTATIONS 

A. GENERAL 

The general design of this section is to utilize all of 
the computational concepts discussed in Sections III, IV and 
V in the computations for a single time-step in the heating 
model developed for use in the FNWC prediction model. After 
testing Fo” and ALB with the corresponding values of Raschke 
et al (1973), it was decided that only the computations by 
the 2/3-CL parameterization would be displayed in the meri- 
dional cross-sections (Figs. 6,7,8,9 and 10). The appropriate 
calculations were performed at each gridpoint of the four 
meridians used for presentation pruposes in this section. 
However, the same type of computations using the full-CL 
parameterization were also made, but not presented in this 
section, 
B. GEOGRAPHICAL REPRESENTATION OF THE HEAT-BALANCE 

DISTRIBUTION 

The FNWC gridpoint processed analyses for OOOOGMT, 16 
April 1974 were used at the three Pacific cross-sections, 
ante that for 12006MT, 16 April- 1974, vere used for the 
single Atlantic meridian. This vas done so that the set of 
gridpoints vas considered to be subiect to the actual radia- 
tive-transfer calculations involved for these specific 


25055 Figure 6 depicts in symbolic language the key to 
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the computational entries in Fig. 7, 8, 9 and 10. This 
symbolic list presents the computations made at each radia- 
tive sounding gridpoint (1,J) having data in the form of 
Table I(b). The computations proceed from the top of the 
troposphere to the ocean surface and include the sensible 
and latent heat transfers, by the means of the computations 
of Section V. For purposes of climatological data compari- 
son, Figs. 7, 8, 9 and 10 were developed by interpolating 
gridpoint results to integral multiples of 5-degree incre- 
ments. The interpolation routine to this gridpoint spacing 
made use of the Lagrangian cubic interpolation scheme (after 


Spaeth, 1975). 


E io, (ES MES ED 
NA a (221) (2-3) (2-4) 


(6-1) 
+ Q (I-1) (1-2) (1-4) a (1-1) (1-2) (1-3) 
3 (3-1) (3-2) (3-4) 4 (4-1) (4-2) (4-3) 
Finally for ease in reconciling the magnitudes of all heat- 
transfer rates, the time-dependent solar disposition rates 


have been averaged to 24-hourly rates. 


C. EXPLANATION OF SYMBOLIC TERMS 
1. Cross-Section at Level k=2 (Fig. 6) 
The discussion of all insolation parameters discussed 
277775) .m"Sectüd41on IWdealt with the specific time of day 


that corresponded to the hour angle h for the instantaneous 
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time t under consideration, The ala solar insolation 
dealt with is then 

F(S) = SGÈ)"" cos z Gaz) 

m 

In order to avoid reference to specific map times t, the in- 
Ecantaneous solar hour=angles were h = 35, 10, 55, and 35, 
respectively for cross-sections l, 2, 3 and 4, as depicted 
Mer des. 7, 8, 9 and 10. 

QAVE ‘represents the 24-hour average of F(2) and 
appears as the first input symbol in Fig. 6. Its value is 
considered to be more representative climatologically for 
the data day under consideration than F(2). 


QAVE is derived by the formula 


Cos z 


QAVE = F(2) PEE (6-3) 

where 
dos z= [H SingSing + CosòCos$SinH]/t (6-4) 
Н = агсСоз (-ТапфТапбј (6-5) 


Here ô is the solar declination for 16 April 1974, and H is 
the hour angle at local sunset at latitude $. Cos z in (6-3) 
is equal to the 24-hour average cosine of the zenith angle, 
Eq. 4-2. The 24-hour time averaging period yielding QAVE 
gives heating results consistent in magnitude with the terres- 
trial flux divergences, which change only slightly with the 
time of day. The conversion to expected daily averaged 

solar disposition quantities is compatible with the determi- 


nation of a heat budget for the given date (16 April 1974). 
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Other parameters needed for level k=2 are 


_ Cos z 
 QREF - REYF(t) 5 (6-6) 
vhere 
REF(t)= F(2)- A26- A610-(IA10+ IS10) (6-7) 


REF(t) is the instantaneous solar reflected insolation at a 
gridpoint and QREF is its 24-hour average, assuming that 
planetary albedo remains constant for the 24-hour period, 
although it was calculated for the indicated solar time t. 
This assumption requires that the cloud amounts computed at 
the indicated synoptic times are representative of the 
entire day. 

The same principle will be used with regard to all 
other solar parameters in the conversion from time-dependent 
values at solar time t to 24-hour averaged values. Superior 
bars Gi) are not used in presenting the averaged values 
shown in the cross-sections key, Fig. 6. The 24-hour aver- 


age system balance, BALT, is computed from 
BALT = QAVE - (QREF + F,*) (6-8) 


for the tropopause level k=2 at the indicated latitude. Net 
terrestrial fluxes, such as EF, were considered to be con- 


stant throughout the 24-hour period, a valid assumption if 


the cloud cover remains quasi-constant for the period. 
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2. Cross-Section in Layer (2,6) 

The following definitions apply for the layer (2,6) 
as referred to in Fig. 6. All the heat transfers shown in 
this layer are assumed to be of radiative character only, as 
in the corresponding version of the FNWC heating model. The 


averaged radiative heating (cooling) rate is given by 
BAL26 = Q26 - F26 (6-9) 


where 


Q26 = daily solar absorption in layer (2,6) 


and is defined relative to A26(t) by a cosine transformation 


similar to Eq. 6-3, and 


F26 = terrestrial cooling rate. 
3. Cross-Section in Layer (6,10) 
The 24-hour average radiative cooling in the layer 


(6,10) is given by 
BAL610 = Q610 – F610 (6-10) 


where Q610 is the sum of Q68 and Q810, BAL 68 has been taken 
as one-half of BAL610. However, the heat balance in (8,10) 


has been modified for inclusion of Hr» as given by Eq. 5-7. 


BAL810 = Q810 - F810 + H, = BAL68 + H (6-11) 


Г r 


Here H, is the sensible heat flux convergence in the layer 


r 
2177. .0 15 the süm of F68 and F810, the terrestrial 


flux losses for the respective layers, 
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4. Cross-Section at Air-Sea Interface (k=10) 


The heat balance at the earth's surface (BALB) is as 


defined in the following equation 


= = x 
BALB QABG Fig (Hr + E) (6-12) 


QAGB is the 24-hour average insolation absorbed by the sur- 


face as 
QABG = QABG(t)(Cos z/Cos z) (6-13) 


Hp and E were computed respectively from Eqs. 5-7 and 5-13 
and considered constant for the computational day, 16 April 
1974. 
D. LATITUDINAL CROSS-SECTIONS OF THE VERTICAL HEAT 

BALANCE FOR 16 APRIL 1974 

Figs. 7, 8, 9 and 10 as previously explained represent 
the single time step of heating computations for each of the 
four meridians used in this study. The four figures have 
been divided into (a) tropical results and (b) mid-to-high- 
latitude results. While the results depicted in these cross- 
sections are exhibited as representing daily-averaged values, 
they are actually based upon heat-computations at the speci- 
fic map times of O000GMT and 1200GMT on 16 April 19/4. There- 
fore, for these results to be meaningful as a stepwise part 
of the heat package subroutine of FNWC, the solar radiative 


absorption and reflectance terms would have to be 
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recoverable as a function of GMT, i.e., 


F(2,t) = QAVE* (Cos z/Cos z) (6-14) 
REF(t) = QREF* (Cos z/Cos z) (6-15) 


etc. Thus solar disposition terms may then be utilized in 
connection with the one-hour stepwise application of the 
thermodynamic equation of the set of primitive equations 
used in the FNWC prediction process, assuming the 2/3-CL 
parameterization, 

The full-CL heat-balance cross-sections were also com- 
puted, but have not been presented here for the sake of 
brevity. However, certain key radiative heating (cooling) 
rates have been summarized in zonally-averaged form from 
both cloud parameterizations for the purposes of presenting 
the mean latitudinal comparisons (Section VII) of the 
purely radiative heat budget of the FNWC data as it was 


applicable to 16 April 1974. 
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145E Longitudinal cross-section, tropical section. 
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Figure 9(a). 


Values computed from data for 16 April 1974. 
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145E Longitudinal cross-section, higher latitude section. 
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Figure 101a). 


Values computed from data for 16 April 1974. 
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71 Lobe LATITUDINAL DESTRIBUTION OF 
RADIATIONAL BALANCE TERMS OF THE 
OCEAN-ATMOSPHERE SYSTEM 
A. GENERAL 

The latitudinally-distributed results of the purely 
radiative contributions, over the ocean-atmosphere system, 
are presented in Figs. 12(a,b) for the 2/3-CL and 13(a,b) 
for the full-CL cases respectively. These cross-sections 
show the results after averaging over: the four meridians 
considered in this study. The results are displayed in the 
format of Fig. ll. Radiative transfer Q68 and Q810 have 
been combined with F68 and F810 of Fig. 6 to compute layer- 
average radiational warming rates in the layer (6,10). The 
turbulent transfer heating rates E + Hp into the atmosphere 
(and the corresponding cooling rates in the ocean) which 
were included in the results of Figs. 7, 8, 9 and 10 are now 
omitted in the radiative composites presented in Figs. 12 
amd 13, 

In obtaining latitudinally-distributed means of radia- 
tive heating rates at each five-degree multiple of latitude 
6, all values of each radiative parameter listed in Fig. 11 
at latitude 6 in the range 205,...,65N (5-degree increments) 
were simply averaged over the four meridians providing the 
averaged value at $. At $ = 65N there was only one contri- 


bution to the latitudinal average, while in the Southern 
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Hemisphere latitudes, 205, 15 and 10S, only two values (on 
À = 125W and A = 35W) of each parameter contributed to the 
means listed in Figs. 12 and 13. At 5S, there were three 
sets of radiative parameters entering into Q (4). Otherwise, 
there were four values of each radiative parameter entering 
into the computed radiative mean values Q(6) Of the cross- 
sections, Therefore, near the northern and southern 
boundaries of Figs, 12 and 13, the listed values may not be 
equally representative, although the general trend should be 
reliable. Fig. 14, using data extracted from Figs. 12 and 
13, tends to confirm this statement. 

Figs. 12, 13 and 14 are also useful in determining the 
relative Merits of the 2/3-CL model as compared to the full- 


Mimo de Y. 


B. EARTH-TROPOSPHERE SYSTEM RADIATIONAL BALANCE SUMMARY 

The latitudinally-averaged distributions of the three 
key radiative parameters, Ro» R, and R as functions of lati- 
tude for both the 2/3-CL and full-CL cases are shown inm 
Fig. 14. The radiative-balance parameters were defined 
(after Malkus, 1962) as R (ó), the mean radiative-energy 
transfer rate across the top of the ocean-tropospheric 
System at k=2 (referred to as BALT in Fig. 11); R. (0), the 
mean radiative cooling rate in the troposphere (BAL26 + 
BAL610 in Fig. 11); and R(o), the mean radiative warming 


Cooling) rate at the earth's surface (BAL1O in Fig. ll). 
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2/3-CL model. 
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Similar to Fig. 12(a) except full-CL model (after Smagorinsky, Eq. 2-4) 
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The relationship between these radiative parameters is 
R= REFER (7-1) 


Tables X, XI and XII list the zonally-averaged values of R 
s 


> 


R and R, as a function of latitude è for the 2/3-CL and the 
full-CL cases. Weighted averages for any Q-value are then 


computed using the following cosine weighting scheme (Eq. 7-2, 


Spaeth, 1975). 


18 fè 
X “r fs su 
4 22 ko 557: 
V dE j = 
ie. avg. (7-2) 
18 
> + 
Фи Сосз Q. 
izi 4 ` 
where Qaşı = value on meridian j at latitude P. and k; is the 


number k=1,...,4 of meridional observations available for the 
meridional averages at Фу. Note that i=1l,...,18 corresponds 
to the 18 latitudes, $; = 20,...,65 latitude degrees. Tables 
X, XI and XII depict comparative values of Ro» R and R, by 
latitudes for both 2/3-CL and full-CL parameterizations for 
то tanuayy, 16 April and 16 July 1974. 

In Fig. 14, the dominant net-flux terms across inter- 
faces are R ($) and R(0), The net flux values obtained from 
the 2/3-CL model (denoted Roi)? exceeded those by the full-CL 
model (R 5) at all latitudes. Both Ro functions had maxima 


at 30N on 16 April 1974. Comparative values of the mean 


150,3 





latitudinal distributions of the radiative balance R(¢) at 
sea level for the 2/3-CL and full-CL cases were also ex- 
tracted from Figs. 12 and 13. They were tabulated in the 
center portion of Tables X, XI and XII for 16 April 1974 and 
plotted also as a function of latitude in Fig. 14. It is to 
be noted that the 2/3-CL case (denoted R.) is greater than 


R, over the entire range of latitude under consideration in 


Fig. l4. The difference R; _ Roi is very nearly equal to 


so that the net cooling in air, R = R 


5::- 
5 o al = 


2 S 2° 


over the entire range of latitude in Fig. 14. 
Im summary, for 16 April 1974, Table X and Fig. 14 show 
that a substantial positive difference in net flux exists at 


the troposphere k=2 for the case of R (for 2/3-CL) over 


1 
the case of R (for full-CL). The zonally-averaged R -dif- 


ference for the range of latitudes in Spring is .0380 cal 


“95 nin İ. Likevise, Table XI and Fig, 14 shov that the net 


flux at the surface for the 2/3-CL exceeds the full-CL by 
-1 


.0351 cal euer , when both sets of values R; and R, have 


been meridionally-averaged. However, the difference Е1-Едә 
ии ii giblè regardless of whether 2/3-CL or full-CL para- 
meterization was employed. 

ANG Ontsisted with the Winter (16 January 1974) studied 
by Spaeth, no change in sign in R = Roo аиләви R, - R, 
occurs for $ > 40N. The reason for the consistent sign o f 


.03 cal ein > which exists at 


the difference R - R 
sl s2 
255500) 11 iatitudes in April is that the "summerlike- 


effect" of an increase in effective insolation with decreased 
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Figure 14. Radiational balance at the tropopause (R.),at 
the ocean surface (R), and in the tropospheric column (R_). 
Solid lines denote computations made with the 2/3-CL (sub- 
scripted "1") cloud model and dashed lines correspond to the 
full-CL model, (subscripted "2"). All computations are for 
16 April 1974 in ly nm) 
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cloud-cover is dominant at all $ < 65N at this time of the 
year (when the declination is 8.5 degrees), whereas there 

is a smaller retention of IR net flux by the troposphere when 
the cloud-cover is reduced to 2/3-CL. These effects apply at 
all latitudes under study here, but undoubtedly, there could 
be a change of sign in both R = R 2 and R; - R, in lati- 


tudes further removed from the subsolar latitude for 16 April 


1974, 


C. CROSS-SEASONAL EFFECTS 

Tables X, XI and XII show the cross-seasonal effects in 
Ros R and R, centered on 16 April 1974. It is clear that 
R increases at each @ in the Northern Hemisphere in the time 
sequence 16 January through 16 July 1974. A similar effect 


applies to R (full-CL case) through the same time sequence. 


2 
Similar results apply to the time sequence for the sizable 
increases of net flux (R) at the earth's surface, both for 
the 2/3-CL and full-CL cases. These cross-seasonal results 
tend to confirm the ee of the radiative model, 

In the Southern Hemisphere latitudes of Tables X, XI and 
XII, there is a systematic trend toward decreasing R rvalues 
in the time sequence 16 January through 16 July 1974. The 
Same trend applies to R-values across the time sequence. 
Considering that this sequence corresponds to the periods of 
summer-fall-winter in the Southern Hemisphere, this cross- 


seasonal trend also appears to afford qualitative support to 


5H777 .1.0tional model tested here, 


EL 





D. COMPARISON OF NET FLUX AT THE TOP OF THE MODEL 
ATMOSPHERE WITH SATELLITE OBSERVATIONS 


Based upon observations of Raschke et al (1973) of 
F2(RAS) and of the albedo ALB(RAS) over the latitude-range 
of this study, the net flux parameter at the top of the 
atmosphere was defined for the NIMBUS III data-period (1-15 


May 1969) by Raschke et al (1973), as 


2.00 


= * 
E o 


(1-ALB(RAS)) - F2(RAS) (7-3) 
Raschke's computed values of the RN for this data period 
used a slightly different solar constant than that employed 
in this study (S = 2.00 ly 2500507 Hence for purposes of 
consistency in comparisons, our average extraterrestrial in- 
solation, QAVE*(2.00/1.92) was used in the definition of RN 
immed. 7-3). 

Values or RN were computed after Raschke at each five- 
degree latitude intersection on the four meridians where 
comparisons with the present model computations were possible. 
For the corresponding radiative model computations of this 


study, the following analog to Eq. 7-3 was adapted. 


2.00 
la 





RNMOD = QAVE* (1-ALBMOD) - FF2 (7-4) 


FF2 represents a parameter analogous to F2(RAS), values 
007 465 vere compared in Section III. E. on a latitudinal 
basis. Likewise, albedo values were compared on a lati- 


tGawan mal basis in Section IV. G. 
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The computational values of RNMOD for 16 April 1974 were 
computed using both the 2/3-CL and the full-CL cloud para- 
meterizations. After obtaining averages across the four 
meridians, the latitudinal comparisons are depicted in Table 
XIII. 

The first three columns of Table XIII show that the lati- 
tudinal distribution of RNMOD is somewhat smaller than the 
values of RN (after Raschke) at all latitudes except 55N. 


The difference 
RN - RNMOD 


is generally positive because of the consistently small albedos 
of the Raschke observations, This difference is considerably 
larger for the full-CL case, when the model albedo was 
recognized as being excessively high, 

A comparison of the mean radiative balance R at the 
tropopause k=2 has been made with that at k=0, which has been 
denoted RNMOD. The results have been tabulated as a function 
of latitude in Table XIII. The R values for the earth- 
troposphere have been extracted from Table X and appear in 


dgölumns four and five of Table XIII. Recall that 


5— :h b EF. (7-5) 


Mat camnparing Eqs. 7-4 and /-5, it is clear that the incoming 


insolation used in the computation of Ro has been reduced by 
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4% of an amount which is approximately 


QAVE — QREF = .3602 ly má a (7-6) 
averaged over all latitudes. However, R 1 exceeds RNMOD(2/3) 


by .017 ly Gat ` ` which must therefore be accounted for by 


the fact 


for the 2/3-CL case. Approximately the same difference 
exists between comparisons for the full-CL cases. 
The IR-flux divergence of (7-7) is due to the presence 


of water-vapor and CO, emitter masses operating in conformity 


2 


with the F,* model (Eq. 3-5). These masses have not been 
subjected to compensating solar insolation absorption in the 


stratosphere, and the 4% solar absorption by oxygen and ozone 


is not sufficient to establish a stratospheric balance. 
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Mid-Seasonal Dates 





Bat. 16 Jan 1974 16 April 1974 16 July 1974 
Rel = a 7. sui :7 
205 . 0681 .0216 -.0229 -,0469 m lli liz  —. 1124 
15 .0552 40052 -.0119 -.0418 -.0817 -.1010 
10 .0488 -.0039 .0002 -.0358 -.0546 -.0772 
5 . 0475 . 0046 .01/1 -.0235 -.0428 -.0616 
0 2015 -.0270 1200 —,0079 .0132 -.0171 
5 -.0050 -.0478 . 0283 . 0243 .0297 -.0062 
10 -.0140 -.0428 20196 “0006 . 0521 0155 
kə -.0489 -.0736 430220 “€.0235 . 0635 s0211 
20 -.0813 -.1010 .0204 -.0368 «0627 . 0136 
25 -.0938 -.1058 .0256 -.0126 . 0825 .0380 
30 -.0998 -.1041 . 0547 . 0293 .0945 .0560 
35 -.1497 -.1551 . 0488 .0301 .1018 .0607 
40 2515777“ .1706 .0540 «0296 1254 .0888 
45 -.1884 -.1835 .0427 0096 0884 -0612 
50 -.2062 -—.1969 -.0055 -.0362 . 1461 SS 
55 -,1798 -.1671 «OO16 -.0146 «1063 «0766 
60 2050:7:-1:75? -—.0428 -,0591 .1063 „1092105 
65N -. 2265 — 2265 = 0039 -.0374 1990 . 1640 
Wives. -.1000 =. 1139 .0227 -.0153 .0523 . 0227 


Table X. Cross-seasonal comparison of zonally-averaged 

xalis Rivas a function of latitude. The 
averaged values for "the 2/3-CL case are denoted with a 
subscripted l and the full-CL case values by a subscripted 
277 5:11 values in ly(min)71). 
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Mid-Seasonal Dates 


Lat. io dra s; 16 April 1974 16 July 1974 
Ер Rə E > Pi Bə 
208 mile 00 0313 0514 1050017056 
15 13607 1125 lots 0681 (0954 0852 
10 o 7500 ` 212909" 0057 Бил 1. 
5 1050017715 a 21875 77 M224 
0 1421007 007 106) 1376 2 a 
5 1100 77 0 pass 1202 14015 2163 
10 0651 01 .1460 .2041 .1726 
15 1227710056 0000 .1470 OO Me 
20 11007 0 oe 1569 155710 
25 10000 os — 2079 1769 55001 ə 
30 2020500 2153 1915 15500” 
35 110177 015 1889 1684 Arie P 
40 550077 00 -1824 ¿1549 .2927 ..2605 
45 2270257 0 1670 ¿1343 DAG O 
50 OO .1256 0952 OS 
55 550400 0500 . 1062 .0850 1101071 
60 27007 “663 .0816 .0603 1125. 53 
65N 2000000630 0614 0206 ANS 
Wt.Avg. OS 1657 (1333 ‚1879 .1349 


Table XI. Cross-seasonal comparison of zonally-averaged 

values of R as a function of latitude Q. The 
averagfed values for the 2/3-CL case are denoted with a sub- 
scripted 1 and the full-CL case values by a subscripted 2. 
711 values in ly (min)”1). 
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Mid-Seasonal Dates 


Lat. 16 Jan 1974 16 April 1974 16 July 1974 
Ы > San ə a > 
205 -.0942 -.0913 -.1042 -.0983 — 10114:—.1674 
15 -.1073 -.1075 -.1134 -.1099 7777:5::5:: 
10 VA SEA —. 1297 = 1315 mə1626“—. 1907 
5 -.1426 -.1492 -.1267 -.1284 —. 1804 “1642 
0 u5 7 —.0O9151  -..455 2.10052 216259 
5 “7 77 x16070) —. 1575.) —, 1671 - 1476“ = Z 6 
10 m7517:7:1 5 — 171777 .1824 = 5218321582 
15 =.1863 1962 -.1708 -.1805 = 19.8 2752221656 
20 3177 € 200) -.18142 - =.1937 x.17:035x 52 
25 Se I6 “1823 =.1895 =. 1756 = 7 
30 m 471 —.1606. =,.1622 1 ии 
35 -.1648 -.1696 -.1401 -.1383 = 17577  — 1825 
40 21O77:517:77/  —.1256” -.1253 - ¿1657 21718 
45 75 71768 =; 1243 —.1247 — 1682 5 1691 
50 = 7, 13 -.1318 =. 1286 =; 274 
55 m500777 1071 —.1046" -.0996 2075247 270272 
60 £7755:51: 7 -—.1246 7 -,1194 => 1123 1045 
65N ISSN — 075 0580 -e0650 ADA 
Wt.Av -.1668 _-.1720 _-.1457 -,1486 =. 1626 5500577 











Table XII. Cross-seasonal comparison of the zonally-averaged 
values of R as a function of latitude $. The 
averaged values for the” 2/3-CL case are denoted with a sub- 
scripted 1 and the full-CL case values by a subcripted 2. 


(All values in ly (min)71), 


11: 








RNMOD Ocean-Troposphere RNRAS 

Lat. 2/3-CL Full-CL Е 1 Ro) 
205 “7777 r mə0z29 -.0469 . 05069 
15 = 0394 0700 =. 0119 -.0418 .0554 
10 mş020:2.0566 . 0002 -—. 0358 0745 
0076 -.0272 . 0171 s. 0235 . 1083 
20093 :—03/7 210553 = 00 79 . 1304 
.0134 -.0392 . 0283 20242 «1545 
10 .0115 -.0448 . 01098 -.0364 21277 
15 0172 “=. 0390 0225 - 0555 1337 
20 2016010 0412 . 0204 575 1272 
25 02982 72220350 0256 20176 . 1150 
30 0467 20213 00547 . 0293 «1039 
35 . 0380 . 0193 . 0488 . 0501 07790 
40 .0429 «0185 .0540 0296 20270 
45 0296720036 ә0427 ә0096 20570 
50 = (0122 —, (0430 =. 0055 =. 0362 0203 
55 = ик 203 TA . 0016 -.0146 —u 
60 = 03540". =.0505 -. 0428 2505921 -.0489 
65N – 0483 -.0770 201009 251574 -.0946 
i Avg. 21107 00026 502217 => .0818 


Table XIII. Comparison of the Ocean-Troposphere Net Radia- 
tion R at k=2, with the model computed Net 
Radiation at the top of the atmosphere RNMOD and with 
climatological Net Radiation at the top of the atmosphere 
from Raschke et al (1973) RNRAS. R and RNMOD have been 
computed using both the 2/3-CL and full-CL parameteriza- 
Broens. `R is from 2/3-CL model and Ro is from the full-CL 


model. = ê 





VT EL, ZONALLY-AVERAGED TROPOSPHERIC 
AND OCEANIC HEAT BUDGETS 
FOR 16 APRIL 1974 

A. THE TROPOSPHERIC HEAT BUDGET 

By averaging over the tropospheric columns for the four 
meridians displayed in Figs. 7, 8, 9 and 10 latitude by 
latitude, the latitudinal distribution of E + Hy results were 
obtained and are presented in Fig. 15. Ro» the tropospheric 
radiative net cooling is listed in Table XII. The zonally- 
averaged E + Hp values computed from Eqs. 5-7 and 5-13 are 
displayed as a function of latitude in Table XIV. The re- 


sulting heating rate of the troposphere may then be expressed 


as a function of è by the right side of Eq. 8-1, 


Der, ə ə. = R, + (E + Hp) è (8-1) 
Fat. =20, m ~IO, — 0.0 5. 10. 19% 20. 
E+Hr 2111757207 / 0061597770 16 .0430 0355 20855 .1758 .0645 
EE .25. 30. 39% 40. 45. 50. 53. 60, 65.N 
E+Hr 2?) 5591707 .1499 .1852 .1966 .1318 .1957 .0642 .1314 


Weighted Average = .1032 
Table XIV. Zonally-averaged values of E+H, as a function 
of latitude. (All values in ly (min)-1). 
A storage heating rate of the column, and Qua is the 
required flux-divergence of heat compatible with the heat 


balance equation for the column, for the observed values of 





F E+Hp and Sou Erçurosoiosdepicts each quantity (R. 


a 3 


(E + Hr) and COLPA + 5.)) as a function of latitude for the 
two cloud-model cases. 

As noted in Section VII. B., Е 1 = R over the range 
of latitudes in this study, with little variation in R, as 
a function of latitude. The distribution of E + Hr as a 
function of latitude is independent of cloud-model. There- 
fore, Fig. 15 shows only one curve for E + Hp e The curves 
of OŒ zı S. are also shown in Fig. 15, as the residual of 
the right sides of Eq. 8-1 for both the 2/3-CL and the full- 


CL cases, respectively. Note that as +5 slightly 


7: 
exceeds (Q + S ), from 10S to 30N as a result of the cor- 
va а 2 


responding excess of R over R for the same latitude 


1 2 


range. Note that both we + 5.) distributions reach maxi- 
mum positive values at @ = 50N, with a general positive 
trend from 35N to 55N. This feature follows from Eq. 8-1, 
since 55 + 5.) must follow the trend of (E + Н.) where 
the latitudinal variation of the latter function peaks more 


sharply, relative to the more or less uniform behavior of 


R (9). 
ae 5° би: 
The cosine-weighted value of a + SA = -.,0425 ly min 
per mə column, Q has not been computed in this study, 


va 


hovever, if the Gm + S. is attributed to mean storage- 
cooling DA of the troposphere alone, the daily storage rate 


corresponds to a temperature-change rate given by 

a (GER À 

n 4 Cyd 8-2 
ет 4.1 AP x 1440 f daya ( ) 


mb 








5 





= 20 


15 
Е*Нг 


10 


(Q va +S3) 2 
I 


l 
l (Qva +Sa)y 
l 


V Ra2 
/ 


7 Rat 
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tre 15. Distribution of Rə EHH, and of 75: as a 
function of latitude for the two cloud-model cases. 


(solid curves reoresent the 2/3-CL case and dashed curves 
represent the full-CL case) 
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with AP mb = 800mb in the troposphere. The resultant cooling 
rate (8-2) over the tropospheric depth of 800mb, considering 


Chat zero lateral tlux-divergence applies, is approximately 


oT _ о = l 
(57) = -,31 C (day) 


2 : 
averaged over the mean cm tropospheric column, 


B. THE LATITUDINALLY-AVERAGED HEAT BUDGET OF THE OCEAN 
The latitudinally-averaged heat-budget of the ocean 
water-mass obeys the following thermodynamic relationship 


for both the 2/3-CL and full-CL cases, 


Q., +S = R + (-(E + H,)) (8-3) 


vo 


Note that Eq. 8-3 is the analog of Eq. 8-1 for the atmosphere, 
with -(E + Hp) representing the average turbulent cooling 
rate for the water-mass by mixing across the sea-air inter- 
face. 

EH Oxu 0 the right side of Eq. 8-3 have been computed 
for each latitude and averaged zonally, with R and E + Hr 
shown in Table XI and XIV,respectively, ə + 5. has been 
computed as the residual of these two terms as a function of 
latitude and the entire set of terms in Eq. 8-3 are graphed 
1577116. 16 £or both cloud cases, 

The primary ocean-mass heating-function is R, the radia- 


tive heating rate at the ocean-air interface, As discussed 


mn Section VIT, Би Е, for all latitudes considered in this 


182:2 





study. 
be greater than (Quo + S.) by 
-(E 4 Hp) q = -(E + H>)" This 


ing rate of the ocean-mass, as 
in the 2/3-CL case than in the 


considered in this study. с 


heat-flux divergence at latitude ф. 


conclusion of Von der Haar and 


Merefore, according Го Ед. 8-3, 5 + 5.) 


1 must 
the same amounts since 

means that the internal heat- 
given in Eq. 8-3, is greater 
full-CL case for all latitudes 
is the required mean oceanic 
Fig. 16 supports the 


Oort (1973) which requires 


increased lateral transport-divergences SE in the oceans, 


based upon the recently observed satellite values of the 


global albedo, 


The values for CQ + S) s from this study 


Seem to corroborate the satellite findings. 


The mean value of (9, + 5.)1 


is 


(Q + 5 


vo 2 


computed between 20S-65N 


SO GE ae Í 


which corresponds to a mean heating rate in the water-mass 


column. 
noticeably larger, 


ponding atmospheric effect. 


and of opposite sign, 


It should be noted that this heating rate is 


than the corres- 
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IX. CONCLUSIONS 


The present work continues the use of the radiational 
two-layer cloud model adapted from Spaeth's mid-winter study 
to FNWC gridpoint data valid over the four meridians for 
16 April 1974. The use of the 2/3-CL parameterization, as 
compared to the full-CL after Smagorinsky (1968), for the 
specification of cloud-layer amounts gave the best verifi- 
cation for the computed global albedo as a function of lati- 
tude when compared with observations from NIMBUS III 
satellite climatology comparable to that of 16 April 1974. 
The use of the 2/3-CL parameterization also gave good veri- 
fication of computed long-wave flux to space when compared 
to the latitudinal distribution of the same element as given 
by satellite climatology. 

The chief shortcoming of the computed global albedo 
versus the satellite measured albedo occurred in tropical 
latitudes (20S-20N). Here, it appears that cloud-elements 
do not behave as reflecting sheets for insolation but rather 
as focusing devices through the intra-cloud spaces. The use 
of 50-60% model cloud reflectances seems excessive in these 

/ 
latitudes. As a practical alternative it may be reasonable 
to tune the cloud reflectances until agreement of the com- 
puted global albedos with satellite measurements is 


achieved in the latitude range 20S-2ON, 
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In the topics, all analyses are subject to data-sparsity 
and hence are largely dependent upon climatology and a rela- 
tively few radiosonde reporting stations. An apparent weak- 
ness in this study is that of overspecifying cloud-amounts 
through the use of Smagorinsky type formulas. This could be 
a result of inadequately resolved relative humidities, which 
in general might well be lower than inferred by the FNWC 


analyses, 
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